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Abstract 


This  paper  discusses  an  analytical  performance 
prediction  methodology  for  annular  propulsive  noz¬ 
zles  and  the  application  of  that  methodology  to 
selected  nozzles.  Thrust  efficiencies  and  static 
pressure  profiles  from  cold  flow  testing  of  the 
selected  nozzles  are  sumnanzed  and  compared  to  the 
predictions.  These  comparisons  show  that  the  ana¬ 
lytical  methods  are  reasonably  accurate  if  the 
radial  velocity  components  in  the  transonic  region 
of  annular  nozzles  are  small.  The  areas  of  pre¬ 
diction/test  disagreement  are  identified  for  future 
improvements. 
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jntroQuet ion 

The  conventional  propulsive  nozzle  it  an  j*!- 
symmetric  converging-diverging  design  and  has  been 
used  successfully  in  rocket,  turbojet  and  ramjet 
engines  for  many  years.  'he  advantages  of  sue n 
nozzles  include  efficient  design. point  performance, 
reasonable  weight  and  cost  and  the  availability  c i 
accurate  design/analysis  techniques.  On  the  other 
hand,  the  disadvantages  of  such  nozzles  include 
excessive  length  and  inefficient  oft-desiqn  per¬ 
formance  . 

Several  annular  nozzle  concepts,  formed  tv 
inserting  axisymmetric  centerbodies  within  conven¬ 
tional  nozzles,  have  been  considered  to  overcome 
these  disadvantages.  As  depicted  in  Figure  i,  the 
centerbody  may  be  located  either  upstream,  down¬ 
stream  or  in  the  throat  plane  of  the  conventional 
nozzle.  The  objective  of  the  present  work  was  to 
conduct  analytical  and  experimental  investigations 
of  selected  annular  nozzles  to  acquire  a  better 
understanding  of  the  flow  processes  and  performance 
potential  of  these  concepts. 

This  paper  presents  the  results  of  this  worn 
for  five  annular  nozzles  which  encompass  the  desiqn 
options  noted  in  Figure  1.  Primary  emphasis  is  on 
an  analytical  performance  prediction  methodology 
for  such  nozzles  and  the  application  of  tfat  meth¬ 
odology.  The  results  of  considerate  cold  flow 
testing  of  selected  nozzle  models  are  '  ..marized  and 
compared  with  the  analytical  predict  .-ns. 

All  analytical  performance  '  redictions  pre¬ 
sented  in  this  paper  were  made  r  1 jr  to  nozzle  model 
cold  flow  testing.  No  revise^  predictions  based  on 
post-test,  "fine  tuned'1  reruns  of  the  analytical 
codes  are  included.  Therefore,  the  comparisons  with 
experimental  data  show  both  the  strengths  and  wedK  - 
nesses  of  the  basic  prediction  methodology.  Some 
additional  capabilities  which  would  enhance  the 
prediction  methodology  are  identified  as  areas  for 
future  improvemer  .  $ . 


Cold  Mow  Testing 

AnPw'dr  nozzle  model  cold  flow  testing  was 
conduc  ed  in  the  McOonnell  Aircraft  Company  Nozzle 
Test  Stand  (NTS).  The  nozzle  dr»vpn  ejector  as* 
seep  j  y  was  used  during  testing  to  obtain  performance 
data  at  higher  pressure  ratios.  A  brief  description 
o  the  NTS  and  ejector  assembly  are  presented  in 
this  section. 
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Fig.  1  Annular  nozzle  concepts. 

The  NTS,  depicted  in  Figure  2,  was  developed  to 
provide  precise  measurements  of  nozzle  thrust  and 
flow  rate.  Test  nozzles  are  attached  to  a  plenum 
assembly  which  is  mounted  between  the  arms  of  a 
double  armed  pendulum.  Pressure  regulated,  dry  air 
is  supplied  to  the  plenum  by  dual  supply  1 ines  which 
run  parallel  to  the  pendulum  arms.  Thus,  the 
incoming  air  is  perpendicular  to  the  nozzle  thrust 
axis  and  its  momentum  does  not  influence  the  thrust 
measurement. 

An  ejector  assembly,  also  depicted  in  Figure  2, 
is  used  with  the  NTS  to  increase  nozzle  pressure 
ratio  (plenum  pressure/back  pressure)  capability. 
The  exhaust  flow  from  the  test  nozzle  drives  the 
ejector  and  lowers  the  nozzle  back  pressure  below 
the  ambient  pressure. 

Nozzle  model  thrust  during  testing  is  measured 
by  recording  the  deflection  of  the  pendulum  about 
the  I-beam  flexure  with  a  linear  variable  differ¬ 
ential  transformer  (IVDT).  The  thrust  magnitude  is 
determined  by  multiplying  the  deflection  distance 
by  the  deflection  constant  which  is  determined  prior 
to  each  test  series  using  dead  weight  loading  with  a 
pulley  arrangement.  The  air  supply  lines  are 
pressurized  during  the  calibration  process.  Thus, 
tare  forces  due  to  thrust  and  pressure  loads  are 
included  in  the  calibration  data. 


PLENUM  GAP 
SEAL 

Fig.  2  Nozzle  test  stand. 

The  thrust  measurement  accuracy  of  the  NTS  for 
nozzles  exhausting  into  the  atmosphere  (i.e..  with¬ 
out  the  ejector)  were  verified  by  tests  with  an 
American  Society  of  Mechanical  Engineers  (ASME) 
standard  nozzle.  These  nozzles  are  commonly  used  to 
verify  thrust  measurement  facilities  since  the 
thrust  efficiency  is  easily  calculated  with  one- 
dimensional  theory.  As  shown  in  Figure  3,  a  curve 
fit  of  data  from  a  typical  verification  test  series 
with  the  ASME  nozzle  is  approximately  0.4  percent 
below  the  theoretical  line.  This  excellent  corre¬ 
lation  confirms  the  accuracy  of  the  NTS  thrust 
measurement  method. 
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Fig.  3  NTS  thrust  measurement  evaluation 
surroiary. 

Thrust  measurement  capability  of  the  NTS  with 
the  ejector  was  verified  by  comparing  nozzle  effi¬ 
ciency  data  between  runs  exhaustina  into  the  atmo¬ 
sphere  and  runs  exhausting  into  the  ejector.  Such 
comparisons  for  a  conical  nozzle  are  presented  in 
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Figure  4  and  show  excellent  agreement  Detween  the 
ejector  and  atmospheric  exhaust  runs,  verifying  the 
thrust  measurement  capability  of  the  NTS  with  the 
ejector. 
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Fig.  4  NTS  thrust  measurement  evaluation 
summary. 
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The  propulsive  nozzle  is  the  thrust  producing 
component  of  jet  propulsion  engines  su„n  as  rennets, 
turbojets  and  ramjets.  The  function  o'  the  propul¬ 
sive  nozzle  is  tn  convert  the  random  tnermai  energy 
of  the  gases  entering  the  nozzle  into  directed 
iinetic  energy  at  the  nozzle  exit,  thus  creating  a 
thrust  on  the  propulsive  system.  'he  expansion 
process  must  be  accomplished  as  efficiently  as 
possible,  since  the  gross  thrust  developed  by  the 
nozzle  is  directly  related  to  tne  nozzle  exhaust 
velocity.  For  rocket  engines,  the  gross  thrust 
produced  by  the  nozzle  is  aiso  the  net  thrust. 
However,  for  airbreathing  engines  (turbines  ar.o 
ramjets),  the  net  thrust  is  the  difference  between 
the  nozzle  gross  thrust  and  the  inlet  ram  draa. 
Consequently,  all  losses  in  nozzle  gross  thrust  for 
atrbreathers  are  magnified  several  times  in  net 
thrust  and  it  is  imperative  that  the  nozzie  effi¬ 
ciency  be  as  large  as  possible. 

Figure  6  illustrates  the  thrust  proouemo  com¬ 
ponents  of  a  generic  annular  nozzle,  'he  nozzle- 
thrust  is  composed  of  the  thrust  developed  across 
the  throat  Ty,  the  thrust  developed  bv  the  pressure 
acting  on  the  cowl  surface  If,  the  thrust  developed 
by  the  pressure  acting  on  the  centerbody  surface  T p 
and  the  thrusts  acting  on  the  centerbody  base  area 
Tpg  and  the  annular  base  area  around  the  nozzie  e«ft 
tCB- 


The  mass  flow  rate  through  the  test  nozzle  is 
metered  through  a  choked  venturi  which  was  cali¬ 
brated  by  the  Colorado  Engineering  Experimental 
Station.  Air  temperature  and  pressure  measurements 
in  the  venturi  are  made  to  correct  for  real  gas 
effects  on  mass  flow  rate. 

Four  annular  nozzles  which  encompass  the  de¬ 
sign  options  noted  in  Figure  1  were  selected  for 
model  fabrication  and  cold  flow  testing.  The 
internal  contours  of  these  nozzles  were  developed  by 
various  design  techniques^.  An  additional  nozzle 
model  was  also  fabricated  based  on  the  design 
improvement  study  described  later  in  this  paper. 

Figure  5  summarizes  the  nozzle  model  cold  flow 
performance  testing  conducted  in  the  NTS  facility. 
The  nozzle  model  and  the  number  of  test  runs 
conducted  with  atmospheric  and  ejector  exhaust 
conditions  are  sunmarized.  Twenty  test  runs  were 
made  and  data  were  recorded  at  8  to  10  values  of 
nozzle  pressure  ratio  for  each  run. 
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NOZZLE 

MOOEL 

NUMBER  OF  RUNS 

EXHAUST 

TOTAL 

ATMOSPHERE 

EJECTOR 

N» 

- 

2 

2 

N? 

3 

3 

6 

N3 

7 

3 

5 

N4 

i 

7 

3 

NS 

4 

4 

total  number  of  runs 

20 

Fig.  5  Cold  flow  test  sumary. 
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Fig.  6  Thrust  components  of  an  annular  nozzle. 

The  relative  magnitudes  of  these  thrust  compo¬ 
nents  depend  on  the  nozzle  design.  For  nozzle 
geometric  area  ratios  between  3:1  to  5:1,  as  con¬ 
sidered  in  this  study,  the  magnitude  of  the  throat 
thrust  is  approx imatel y  80  percent  of  the  nozzle 
ideal  thrust,  the  cowl  thrust  is  from  iO  to  15 
percent  of  the  ideal  thrust  and  the  centerbody  and 
base  thrusts  and  losses  comprise  the  remainder. 
However,  if  the  throat  flowfield  has  a  laroe  radi¬ 
ally  inward  or  outward  component,  the  axial  compo¬ 
nent  of  the  throat  thrust  is  reduced  considerably 
and  the  cowl  and  centerbody  thrust  components  become 
a  larger  portion  of  the  total  thrust. 

The  accuracy  and  efficiency  of  nozzle  perfor¬ 
mance  prediction  techniques  are  highly  dependent  on 
the  flowfield  model  chosen  to  represent  the  actual 
flowfield.  In  the  present  study,  the  flowfield 
model  is  based  on  the  foi lowing  assumptions: 

(1)  steady  axi symmetric  flow. 

(2)  invisetd  nonconductina  fluid  with  no  body 
forces . 

(3)  thermally  and  calorically  perfect  qas  and 

(4)  separation  and  base  pressures  selected 
from  empirical  correlations. 
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The  governing  equations  consist  of  the  continuity 
equation,  the  component  momentum  equations,  the 
energy  equation  and  the  thermal  and  caloric  equa¬ 
tions  of  state.  A  detailed  discussion  of  these 
equations  is  presented  in  Reference  3. 

The  choice  of  a  numerical  method,  or  methods, 
to  solve  these  governing  equations  is  also  crucial 
to  the  success  of  the  performance  prediction  pro¬ 
cedure.  A  suDsonic/transomc  solution  is  required 
to  define  the  throat  thrust,  while  a  supersonic 
solution  is  required  to  define  the  cowl,  centerpody 
and  base  thrusts.  Experience  has  shown  that  sub- 
somc/transonic  flowfields  can  be  adequately  deter¬ 
mined  by  finite  difference  approaches,  whereas  the 
method  of  characteristics  is  the  most  accurate 
approach  to  determining  supersonic  flowfields  in 
nozzles.  Therefore,  the  numerical  methods  selected 
for  the  performance  prediction  procedure  in  the 
present  study  are  a  finite  difference  approach,  as 
employed  in  the  VNAP  code4,  to  define  the  sub¬ 
sonic/transonic  flowfield  and  a  method  of  charac¬ 
teristics  approach,  denoted  as  the  MOC  code^,  to 
define  the  supersonic  flowfield. 


Performance  Prediction  Example 

The  application  of  the  se'ected  performance 
prediction  methodology  to  nozzle  Nl,  see  Figure  7, 
is  presented  in  this  section.  The  flowfield  solu¬ 
tions  obtained  from  the  VNAP  and  MOC  codes  are 
presented  and  the  calculations  of  throat,  cowl, 
centerbody  and  base  thrusts  are  described  in  detai  1 . 
Comparisons  with  cold  flow  test  data  are  also  made. 
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Fig.  7  Geometry  of  nozzle  Nl. 

Figure  7  depicts  the  geometry  of  nozzle  Nl. 
The  nozzle  has  an  axisymmetric  cowl  consisting  of  a 
cylindrical  upstream  section,  a  conical-circular 
arc  shoulder  section  and  a  contoured  downstream 
section.  The  axisymmetric  centerbody  has  a  cylin¬ 
drical  upstream  support  section,  a  conical  initial 
section,  a  circular  arc  shoulder,  a  conical  section 
downstream  of  the  shoulder  and  a  flat  aft  facinq 
base.  Note  that  the  shoulder  sections  of  the  cowl 
and  centerbody  are  at  the  same  axial  station. 

As  previously  mentioned,  the  subsonic/tran¬ 
sonic  solution  obtained  from  the  VNAP  code  defines 


the  nozzle  throat  flowfield.  Figure  8  illustrates 
the  Mach  number  distribution  in  the  throat  region  of 
nozz le  Nl  as  an  ex  amp (e  of  t he  f 1 owf ie 1 d  definition 
provided  by  the  VNAP  code,  from  this  definition, 
the  throat  thrust  is  determined  by  integrating  tne 
momentum  flux  and  pressure  forces  across  a  selected 
line  that  spans  the  throat  in  the  slightly  super¬ 
sonic  region.  This  same  line,  denoted  as  the 
initial  value  line  (IVL),  is  then  used  to  start  the 
MOC  code  solution  of  the  supersonic  flowfield. 
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Fig.  8  Mach  number  distribution  for  nozzle  Nl. 

The  thrusts  developed  by  the  pressures  acting 
on  the  cowl  and  centerbody  surfaces  are  determined 
by  integrating  the  axial  components  of  the  pressure 
forces  acting  on  these  surfaces.  The  pressure 
distributions  along  the  cowl  and  centerbody  sur¬ 
faces  are  obtained  from  the  supersonic  flowfield 
solution  provided  by  the  MOC  code.  The  MOC  code 
constructs  the  flowfield  solution  along  the  pre¬ 
viously  mentioned  IVE,  beginning  at  the  cowl  and 
proceeding  along  the  IVL  to  the  centerbody  surface. 
The  code  continues  the  solution  along  the  centerbody 
surface  until  a  specified  separation  pressure  is 
reached.  At  this  point,  a  tree-pressure  boundary  is 
calculated  and  projected  to  the  nozzle  axis.  Left¬ 
running  Mach  lines  emanating  from  these  boundaries 
are  propagated  across  the  flowfield  to  the  cowl. 
The  Mach  lines  are  allowed  to  cross  and  fold  over, 
simulating  the  occurrence  of  embedded  obliaue  shock 
waves.  This  procedure  is  continued  unti I  the  entire 
cowl  and  centerbody  wall  pressure  distributions 
have  been  determined.  Figure  9  illustrates  the 
left-running  Mach  lines  in  the  supersonic  flowfield 
of  nozzle  Nl  as  an  example  of  the  MOC  code  solution. 
Note  the  occurrence  of  a  mild  riqht-runmng  oblique 
shock  wave  near  the  end  of  the  cowl  in  this  nozzle. 

An  alternative  to  the  throat  flowfield  defined 
by  the  VNAP  code  is  to  simply  assume  the  qeometrv  of 
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Fig.  9  Left-running  Mach  lines  for  nozzle  Nl. 

a  slightly  supersonic,  constant  Mach  number  line 
that  spans  the  nozzle  throat.  This  uniform  IVL  is 
then  treated  in  the  same  manner  as  the  VNAP-defined 
IVL  for  throat  thrust  calculations  and  MOC  code 
usage.  As  will  be  shown,  this  alternative  procedure 
provides  good  results  if  the  throat  ''lowfield  is 
reasonably  uniform. 

The  predicted  and  measured  cowl  pressure  dis¬ 
tributions  are  presented  in  Figure  10.  The  pre¬ 
dicted  pressures  are  presented  for  two  I  Vi's:  the 
IVL  determined  from  the  VNAP  code  and  a  uniform  IVL 
with  a  constant  Mach  number  of  1.02.  The  two 
predicted  pressure  distributions  are  reasonably 
close,  indicating  that  the  subsomc/transonic  flow- 
field  for  this  nozzle  does  not  differ  greatly  from  a 
uniform  flowfield.  The  Mach  number  distributions 
previously  presented  in  Figure  8  support  that  con¬ 
clusion.  The  measured  pressure  distribution  agrees 
well  with  the  predicted  results  except  in  the  region 
just  downstream  of  the  cowl  shoulder  where  the 
predicted  pressure  rise  is  more  rapid  than  the 
measured  rise. 

The  cowl  and  centerbody  base  thrusts  are  com¬ 
puted  by  assuming  that  a  uniform  pressure  acts  over 
these  base  areas.  Base  pressures  must  be  estimated 
from  empirical  correlations.  In  the  present  inves¬ 
tigation,  the  cowl  base  pressure  was  equated  to  the 
ambient  pressure.  Based  on  a  limited  amount  of 
experimental  data,  the  centerbody  base  pressure  was 
assumed  to  be  five  percent  of  the  nozzle  chamber 
stagnation  pressure. 

The  predicted  and  measured  values  of  the  nozzle 
pressure  ratio  efficiency  are  presented  in  Fig¬ 
ure  11.  (Nozzle  efficiency  definitions  are  dis¬ 
cussed  in  the  Appendix  of  this  paper.)  The  observed 
trend  of  measured  efficiency  as  a  function  of  nozzle 
pressure  ratio  fol lows  the  predicted  trend  very 
closely.  The  magnitude  of  the  predicted  efficiency 
is  approximately  one  percent  less  than  the  measured 
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Fig.  10  Cowl  pressures  for  nozzle  Nl. 
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Fig.  11  Efficiency  comparison  for  nozzle  Nl . 
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efficiency,  which  is  excellent  agreement.  These 
results  demonstrate  that  the  performance  prediction 
metnodology  oeveloped  during  this  investigation  is 
capable  of  accurately  predicting  the  flowfield  and 
performance  of  annular  nozzles. 


Resul ts 

The  results  of  the  flowfield  analysis,  perfor¬ 
mance  prediction  and  experimental  testing  of  three 
other  annular  nozzle  configurations  are  presented 
in  this  section.  These  three  nozzles,  which  have 
their  centerbodies  located  downstream,  upstream  and 
in  the  throat  plane,  exercise  the  capabilities  of 
the  methodology  and  illustrate  its  strengths  and 
weaknesses. 

Nozzle  N2 

Figure  12  depicts  the  geometry  of  nozzle  N2. 
The  nozzle  has  an  ax  1  symmetric  cowl  consisting  of  a 
cylindrical  upstream  section,  a  conical-circular 
arc  shoulder  section  and  a  contoured  downstream 
section.  The  ax i symmetric  centerbody  has  a  cylin¬ 
drical  upstream  support  section,  a  conical  initial 
section,  a  circular  arc  shoulder  section,  a  conical 
section  downstream  of  the  shoulder  and  a  flat  aft 
facing  base.  Note  that  the  centerbody  shoulder  is 
downstream  of  the  cowl  shoulder. 
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Figure  13  presents  the  Mach  number  distribu¬ 
tion  obtained  by  the  VNAP  code  for  this  nozzle.  The 
flowfield  is  highly  distorted  in  the  transonic 
region.  The  Mach  number  contours  show  a  rapid 
expansion  around  the  cowl  shoulder,  with  the  Mach 
number  reaching  values  larger  than  1.5.  This  highly 
supersonic  region  is  followed  by  a  rapid  compression 
which  causes  the  flow  to  go  subsonic  on  the  cowl. 
This  compression  is  caused  by  the  rapid  inward 
turning  of  the  cowl  contour,  coupled  with  the 
considerable  radially  outward  momentum  in  the  tran¬ 
sonic  flowfield. 

The  VNAP  code  solves  the  Nav ler-Stokes  equa¬ 
tions  for  unsteady  flow  using  MacCormack's  finite 
difference  scheme  and  the  steady-state  flowfield  is 
obtained  as  the  asymptotic  solution  at  large  time. 
The  VNAP  code  uses  this  scheme  to  obtain  the 
flowfield  solution  on  lines  which  are  perpendicular 
to  the  axis  of  the  nozzle.  This  approach  works  best 
if  the  finite  difference  grid  is  well  matched  to  the 
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Fig.  13  Mach  number  distribution  for  nozzle  N2. 

physical  geometry  to  be  analyzed.  That  is.  accurate 
flowfield  solutions  are  obtained  if  the  flow  stream¬ 
lines  are  nearly  perpendicular  to  the  solution 
lines.  However,  if  large  radial  flows  are  present, 
the  flow  streamlines  cross  the  solution  lines  at 
small  angles  and  large  errors  in  flow  properties  can 
occur.  Therefore,  the  throat  thrust  calculated  for 
the  VNAP-defined  IVL  for  nozzle  N2  may  be  in 
considerable  error. 

When  the  VNAP-defined  IVL  for  nozzle  N2  is 
employed  in  the  MOC  code  to  determine  the  supersonic 
flowfield,  the  flow  goes  subsonic  on  the  cowl  and 
the  code  aborts  at  an  axial  location  corresponding 
to  the  subsonic  pocket. 

Attempts  were  made  to  obtain  a  supersonic 
solution  by  starting  the  MOC  code  with  a  uniform  IVL 
across  the  nozzle  geometric  throat  {i.e.,  the  conic 
section  defined  by  the  line  connecting  the  center  of 
curvature  of  the  cowl  shoulder  and  the  center  of 
curvature  of  the  centerbody  shoulder)  for  Mach 
numbers  of  1.01  and  1.02.  Both  attempts  resulted  in 
MOC  code  aborts  due  to  the  computed  flow  becoming 
subsonic  on  the  cowl.  However,  for  a  Mach  number  of 
1.03,  the  flow  remained  supersonic  and  a  flowfield 
calculation  was  obtained. 

Figure  14  presents  the  pressure  distributions 
on  the  cowl  obtained  from  the  MOC  solution  for  the 
uniform  IVL  and  the  VNAP  IVL.  For  the  uniform  IVL 
solution,  the  initial  expansion  is  followed  imme¬ 
diately  by  a  strong  compression  which  decelerates 
the  Mow  to  only  a  slightly  supersonic  Mach  number. 
Downstream  of  the  nearly  constant  pressure  reqion, 
the  pressure  drops  smoothly.  For  the  VNAP  IVL,  the 
flow  initially  expands  to  a  lower  pressure  because 
of  the  effects  of  the  upstream  qeometrv.  The 
recampresswn  then  results  in  a  subsonic  flow  region 
that  causes  the  MOC  code  to  abort.  The  experimental 
measurements  of  static  pressure  on  the  cowl,  also 
presented  in  Figure  14,  indicate  that  the  initial 
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Fig.  14  Cowl  pressures  for  nozzle  N2. 

expansion  is  more  rapid  than  the  uniform  I VL  solu¬ 
tion  predicts,  but  not  as  rapid  as  the  VNAP  IVl 
solution  indicates.  The  recompression  is  clearly 
indicated,  but  the  flow  apparently  remains  super¬ 
sonic.  Note  that  the  computed  and  measured  pres¬ 
sures  agree  well  at  more  downstream  cowl  locations. 

The  predicted  and  measured  nozzle  efficiencies 
for  nozzle  N2  are  compared  in  Figure  15.  The  dual 
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Fig,  15  Efficiency  comparison  tor  nozzle  N2. 


values  of  measured  efficiencies  at  lower  pressure 
ratios  are  due  to  centercody  wake  structure  cnanges 
wnich  were  determined  by  exit  plane  pitot  pressure 
surveys  not  discussed  in  this  paper.  At  mgner 
pressure  ratios,  the  predicted  efficiencies  o'* 
approximately  2  percent  higher  than  me  measured 
values.  This  difference  could  be  anticipated  since 
the  predicted  cowl  pressures  (Figure  14!  are  gen¬ 
erally  higher  than  the  measured  pressures.  The 
inability  to  more  accurately  analyze  the  rlowfieid 
is  disappointing,  especially  since  the  predicted 
presence  of  the  high  pressure  region  suggests  that 
this  type  of  nozzle  could  be  modified  to  a  higher 
performing  design.  This  inadequacy  in  predictive 
capability  is  identified  as  a  major  deficiency  in 
the  state-of-the-art  of  nozzle  1 1 ow field  analysis. 

Nozzle  N3 


figure  16  depicts  the  geometry  of  nozzle  N3. 
The  nozzle  has  an  ax.isymmetric  cowl  consisting  or  a 
cylindrical  upstream  section,  a  con ica 1 -c i rcu 1 ar 
arc  shoulder  section  ana  a  conical  downstream  sec¬ 
tion.  The  ax lsymnetric  centerbody  has  a  cylindrical 
upstream  support  section,  a  conical  initial  sec¬ 
tion,  a  circular  arc  snouider  section,  a  contoured 
section  downstream  of  the  snouider  and  a  f  I  at  alt 
facing  base.  Like  nozzle  Nl.  the  shoulders  of  the 
centerbody  and  cowl  are  at  the  same  axial  station. 
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Fig.  16  Geometry  of  nozzle  N3 

The  Mach  number  contours  and  Mach  lines  ob¬ 
tained  by  the  VNAP  and  MOC  codes  are  not  presented 
since  they  are  very  similar  to  those  previously 
shown  in  Figures  8  and  9  for  nozzle  Nl.  f  iqure  17 
compares  experimental  and  predicted  pressure  dis¬ 
tributions  on  the  cowl  for  both  the  uniform  IVL  and 
the  VNAP  IVL  solutions  with  the  HOC  code.  The 
predicted  pressure  distributions  are  almost  iden¬ 
tical  and  compare  very  well  with  the  measured 
pressures  for  higher  pressure  ratios.  However,  as 
shown  in  figure  17,  the  measured  pressures  at  lower 
pressure  ratios  are  higher  than  the  predicted  val¬ 
ues.  This  difference  is  attributed  to  flow  separa¬ 
tion  from  the  cowl.  This  inability  to  predict  flow 
separation  and  compute  the  resuitimi  ilowfield  is 
another  major  def iciency  in  state-of-the-art  nozzle 
flowfield  analysis. 

Figure  18  compares  the  measured  and  predicted 
nozzle  efficiencies  tor  nozzle  NJ.  uurpr i s inq I v , 
the  predicted  values  are  less  than  the  measured 
efficiencies,  but  the  agreement  is  excellent  at 
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Figure  1.  The  resulting  geometry  for  no22le  N4  is 
illustrated  in  Figure  19. 
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fig.  19  Geometry  of  nozzle  N4. 

The  transonic  Mach  number  distribution  ob¬ 
tained  from  the  VNAP  code  is  presented  in  Figure  ? 0 
and  depicts  a  highly  distorted  flowfield.  A  region 
of  locally  supersonic  flow  occurs  at  the  centerbody 
shoulder.  The  flow  downstream  of  this  reaion 
recompresses  and  remains  suosomc  all  the  way  to  the 
nozzle  axis. 


Fig.  17  Cowl  pressures  for  nozzle  N3. 
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Fig.  18  Efficiency  comparison  for  nozzles  N3 
and  N4. 

higher  nozzle  pressure  ratios.  The  increased  dif¬ 
ference  between  the  measured  and  predicted  effi¬ 
ciencies  at  lower  nozzle  pressure  ratios  merely 
reflects  the  previously  mentioned  flow  separation 
phenomena. 

Nozzle  N4 

Nozzle  N4  is  a  rearrangement  of  nozzle  N3  model 
components  which  locates  the  centerbody  shoulder 
upstream  of  the  cowl  shoulder,  'his  conf igurat ion 
completes  the  annular  nozzle  design  options  noted  in 
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Fig.  30  Mach  number  distribution  for 
nozzle  N4. 

Based  on  simple  one-dimensional  isentropic 
flow  theory^,  it  was  assumed  that  the  flow  would 
choke  at  the  minimum  aeometric  flow  area.  This  area 
is  a  conic  surface  defined  by  a  line  connecting  the 
center  of  curvature  of  the  cowl  shouider  with  the 
edge  of  the  base.  Since  it  was  necessary  to  assume 
an  extension  to  the  centerbody  for  the  VNAP  code 
analysis,  it  was  hypothesized  that  the  assumed 
extension  geomrf'v  influenced  the  results.  Several 
different  center oody  extension  contours  were  eval¬ 
uated  in  an  attempt  to  obtain  supersonic  flow  on  the 
centerbody  or  its  extension.  However,  witi  the 
exception  of  the  previously  mentioned  supersonic 
pocket  at  the  shouider,  the  flow  remained  subsonic 
all  f  he  way  to  the  nozzle  axis  for  reasonable 
extension  geometries. 
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A  slightly  supersonic  !vl  was  selected  from  the 
v'NA?  code  results  end  used  to  start  the  MOC  code 
solution  of  the  supersonic  r  !ow» ield.  'ne  predicted 
pressure  distribution  on  tne  cowl  is  compared  witn 
the  measured  pressures  in  Figure  el.  The  predicted 
pressures  m  the  supersonic  region  downstream  or  the 
cowl  shoulder  do  not  show  the  degree  oi  expansion 
indicated  by  the  experiment  a !  sai-es. 


Fig.  21  Cnwl  pressures  for  nozzle  N4. 

Attempts  were  made  to  start  the  MOC  code 
solution  with  a  uniform  1VI  corresponding  to  the 
previously  mentioned  minimum  flow  area.  These  were 
unsuccessful  due  to  the  computed  ftow  becoming 
subsonic.  As  an  alternative,  the  VNAP  code  was  then 
used  to  obtain  a  solution  for  the  entire  nozzle 
flowfield.  Figure  21  shows  that  the  pressures 
predicted  by  this  VNAP  solution  agree  very  well  with 
the  measured  values  in  the  subsonic  region  upstream 
of  the  cowl  shoulder.  However,  in  the  supersonic 
flow  region  downstream  of  the  shoulder,  the  agree¬ 
ment  is  poorer  than  the  predictions  obtained  with 
the  MOC  code  using  the  VNAP  IVl. 

Figure  18  compares  the  measured  and  predicted 
nozzle  efficiencies  for  nozzle  N4.  The  discon¬ 
tinuity  in  the  measured  efficiency  curve  near  a 
pressure  ratio  of  8:1  is  due  to  flow  separation.  At 
higher  pressure  ratios,  the  predicted  efficiencies 
are  approximately  3  percent  larger  than  measured. 
This  is  attributed  to  the  mismatch  between  the 
computed  and  measured  cowl  pressure  distribution, 
as  previously  discussed.  This  difference  indicates 
that  a  considerable  improvement  could  be  made  in  the 
analysis  caoabi  1  i  t y  for  this  type  of  annular  nozz le. 


Nozzle  Efficiency  Comparison 

Figure  22  compares  the  measured  efficiencies 
of  these  four  annular  nozzles  and  the  reference 
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Fig.  22  Nozzle  efficiency  c •■nrpar i son . 

conical  nozzle.  As  shown,  the  measured  e»  t  »c lent  tes 
of  nozzles  N2  and  fW  are  considerably  lower  tnan 
those  of  tne  conical  nozzie.  However,  the  measured 
efficiencies  of  nozzles  N1  and  N3  indicate  that  they 
could  be  competitive  with  the  conical  nozzle  in  some 
appl ications. 


Oesiqn  improvement  Study 

One  of  the  primary  benefits  of  detailed  anal¬ 
ysis  of  nozzle  flowfields  is  to  identify  means  ter 
performance  improvements.  The  analysis  provides  an 
increased  understanding  of  the  flowfield.  ‘rrm 
which  design  changes  for  improved  performance  mav  be 
proposed.  Analysis  also  provides  an  inexpensive 
alternative  to  experimental  conoarisons  of  proposed 
design  variants.  As  an  example  of  Such  an  analysis 
application,  this  section  summarizes  a  limited 
study  to  improve  the  performance  of  nozzle  N3 
throuqh  small  changes  to  the  suoersomc  contours  of 
the  cowl  and  centerbody.  Nozzie  N3  we-  selected 
since  it  exhibited  good  predicted  performance  and  a 
high  degree  of  confidence  exists  >n  the  analytical 
methods  for  this  configuration.  The  subsonic  con¬ 
tour  of  the  nozzle  was  not  altered  in  this  study  and 
all  computations  were  made  using  the  MOC  code.  A 
uniform  IVL  with  a  Mach  number  of  1.0?  at  the 
minimum  geometric  area  was  used  to  start  the  MOC 
code  solution. 

The  effects  of  three  design  cnanqps  were  in¬ 
vestigated.  They  included: 

(a)  Replacing  the  aft  contoured  centerbody 
surface  with  a  conical  surface, 

(b)  Replacing  the  conical  row!  with  a  Para¬ 
bolic  cowl  and 

(c)  Inserting  a  short  cylindrical  section 
into  the  centerbody  just  downstream  o’ 
the  shoulder. 

The  predicted  results  of  these  desion  chanaes 

were : 

(a)  Replacing  the  aft  contoured  section  o’ 
the  centerbody  with  a  conrji  section 
would  result  in  a  small  performance  de¬ 
crease  . 


9 


(b)  Contouring  the  cowl  boundary  would  result 
in  a  slight  performance  increase  (about 
0.2  percent) . 

(c)  Adding  a  short  cylindrical  section  to  the 
centerbody  would  result  in  a  modest  per¬ 
formance  increase  (about  0.8  percent). 

The  latter  approach  was  selected  for  testing 
since  it  offered  the  greater  potential  for  perfor¬ 
mance  improvement.  This  nozzle  is  denoted  as  NS  and 
is  depicted  in  Figure  23.  The  measured  thrust 
efficiencies  for  nozzle  N5  and  those  for  the  refer¬ 
ence  nozzle  N3  are  compared  in  Figure  24.  This 
comparison  verifies  the  predicted  performance  in¬ 
crease  and  confirms  that  analytical  prediction 
techniques  can  be  used  to  assess  design  changes  for 
performance  improvements. 


CVLINORtCAL  SECTION 
ADDED  AT  SMOULDER 
E  OB  NOZZLE  NS 


Fig.  23  Geometry  of  nozzle  N5. 
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NOZZLE  PRESSURE  RATIO 

Fig.  24  Efficiency  improvement  with 
centerbody  modification. 

The  payoff  in  adding  the  cylindrical  section  to 
the  centerbody  can  be  explained  in  terms  of  cowl 
pressures  by  observing  the  characteristic  wave 


pattern,  aiso  illustrated  in  Figure  23.  fne  pres¬ 
sure  distribution  on  the  cowl  is  primarily  deter¬ 
mined  by  expansion  waves  emanating  irom  tne  center- 
body  surracp  just  downstream  of  Characteristic  Cl. 
which  originates  at  the  intersection  of  the  center- 
body  and  the  I  VC.  These  strong  expansion  waves, 
represented  by  Characteristic  C2  in  Figure  23,  are 
generated  by  the  rapid  turning  of  the  flow  around 
the  centerbody  shoulder  and  cause  a  rapid  decrease 
in  tne  pressure  on  the  cowl.  Adding  a  short 
cylindrical  section  to  the  shoulder  of  the  center- 
body  delays  the  turning  of  the  flow  ana  thus 
translates  the  attendant  pressure  decrease  on  the 
cowl  further  downstream.  Further  study  of  the 
characteristic  pattern  in  Figure  23  also  reveals  the 
limiting  length  of  the  cylindrical  section  to 
achieve  performance  benefits.  Characteristic  C3, 
which  intersects  the  end  of  tbe  cowl,  emanates  from 
the  centerbody  at  a  flow  angle  of  approximately  15 
degrees.  Thus,  changes  to  tne  centerbody  downstream 
of  this  point  have  no  effect  on  the  cowl  pressure. 
Any  other  features  of  the  flow,  specifically  the 
subsonic-transonic  flowfield  and  tbe  pressure  dis¬ 
tribution  on  the  centerbody.  should  not  be  signifi¬ 
cantly  altered  by  this  modification. 


Summary  and  Conclusions 

The  analytical  methods  described  in  this  paper 
are  very  accurate  for  conventional  converging- 
diverging  nozzles  and  are  reasonably  accurate  for 
annular  nozzles  if  the  radial  velocity  components  in 
the  transonic  region  are  small.  The  analytical 
methods  provide  valuable  quantitative  information 
that  can  be  effectively  used  to  understand  the 
flowfield  and  generate  ideas  for  design  improve¬ 
ment.  Furthermore,  the  methods  are  sufficiently 
accurate  that  the  effects  of  small  differences  in 
designs  can  be  accurately  evaluated,  thus  reducing 
the  amount  of  experimental  work  required  to  optimize 
a  given  design  concept. 

Hc-iever ,  there  are  some  deficiencies  in  the 
analytical  capability  for  some  annular  nozzle  flow- 
fields.  These  include  inadequate  transonic  flow 
analysis  capability  for  flows  with  large  radial 
velocity  components,  inability  to  treat  strong 
shock  waves  (particularly  when  the  shock  waves 
produce  subsonic  pockets  in  the  supersonic  flow- 
field),  very  poor  base  pressure  prediction  methods 
and  inability  to  predict  and  calculate  separated 
flows.  Future  improvements  in  these  areas  would 
enhance  the  analytical  capability  for  annular  noz¬ 
zles. 

Appendix 

Nozzle  Efficiency  Definitions 

Due  to  the  importance  of  nozzle  thrust  on 
propulsion  system  performance,  a  consistent  nozzle 
efficiency  definition  must  be  employed  when  compar¬ 
ing  the  performance  of  different  nozzle  designs. 
Since  definitions  of  nozzle  efficiency  in  the  lit¬ 
erature  differ,  the  objective  of  this  appendix  is  to 
present  a  clear  exposition  of  the  various  efficiency 
definitions,  to  careful ly  def me  the  various  param¬ 
eters  that  enter  into  these  definitions  and  to 
recommend  a  consistent  procedure  fcr  the  propulsive 
nozzle  efficiency  definitions  so  that  meaningful 
comparisons  can  be  made  between  various  nozzles. 

Three  different  nozzle  efficiencies  can  be 
def i ned .  They  are : 
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1.  n  [--the  internal  efficiency,  which  compares 
the  actual  internal  performance  of  the  nozzle  (i.e., 
vacuum  performance)  to  the  ideal  internal  perfor¬ 
mance  of  the  nozzle. 

2.  n^R'-the  area  ratio  efficiency,  which  compares 
the  actual  performance  of  the  nozzle  operating  in  an 
atmospheric  pressure  environment  to  the  ideal  per¬ 
formance  of  the  nozzle  operating  in  the  same  envi¬ 
ronment. 


for  DeLaval  nozzles,  the  geometric  iriroai  a tea 
is  defined  as  tne  planar  area  at  the  nozzle  minimum 
cross  section  ana  values  yt  discharge  coefficient 
are  quite  close  to  unitv.  For  ■^conventional 
nozzles,  the  choice  of  geometric  throat  area  is  less 
obvious  and  values  of  discharge  caetricient  may 
differ  consiaeraoly  from  unity.  \n  fact,  tup 
determination  of  the  geometric  minimum  area  for 
annular  nozzies  such  as  discussed  in  tnis  paper  may 
be  quite  difficult. 


3.  npR--the  pressure  ratio  efficiency,  which  com¬ 
pares  the  actual  performance  of  the  nozzle  operating 
in  an  atmospheric  pressure  environment  to  the  ideal 
performance  of  an  adapted  nozzle  with  optimum  ex¬ 
pansion  operating  in  the  same  environment. 

The  effect  of  the  atmospheric  pressure  envi¬ 
ronment  on  performance  is  determined  by  the  nozzie 
pressure  ratio  (NPR),  which  is  defined  as: 

NPR  -  P0/Pa  (A-l) 

where  P0  is  the  nozzle  stagnation  pressure  and  Pa  is 
the  ambient  pressure. 

Each  of  the  three  nozzle  efficiency  defini¬ 
tions  has  some  merit.  The  internal  efficiency 
measures  the  thermodynamic  efficiency  of  the  inter¬ 
nal  expansion  process  and  is  independent  of  the 
nozzle  pressure  ratio.  The  area  ratio  efficiency 
includes  the  thrust  loss  due  to  the  ambient  pressure 
acting  on  the  backside  of  the  nozzle.  The  pressure 
ratio  efficiency  includes  the  loss  in  performance 
when  the  nozzle  is  not  adapted  for  operation  at  a 
particular  nozzle  pressure  ratio. 

The  reference  ideal  thrusts  required  in  these 
efficiency  definitions  also  require  careful  speci¬ 
fication  of  the  reference  ideal  choking  area  and  the 
reference  ideal  exit  area  of  the  nozzle. 

Figure  A-l  illustrates  the  reference  ideal 
choking  area  A*.  It  is  the  area  required  to  pass  the 
actual  mass  flow  rate  and  is  smaller  than  the 
geometric  minimum  (throat)  area  At.  Those  two  areas 
are  related  by  the  nozzle  discharge  coefficient  Cg, 
which  is  defined  as: 


The  actual  mass  flow  rate  may  oe  determineo 
experimentally,  predicted  analytically  or  deter¬ 
mined  from  empirical  discharge  coefficient  correla¬ 
tions.  In  this  study,  reference  mass  flow  rate  was 
determined  using  conventional  steady  one-dimen¬ 
sional  isentropic  nozzie  flow  theory'  ana  actual 
mass  flow  rate  was  predicted  using  tne  Viscous 
Nozzle  Analysis  Program  (VNAP)  code4. 

figure  A-l  also  illustrates  the  reference 
ideal  exit  area  of  the  nozzie.  It  can  oe  defined  in 
several  ways,  for  example: 

(1)  the  nozzle  maximum  internal  area  Api,  or 

(2)  the  nozzle  maximum  internal  area  plus  the 
cowl  base  area  Ae?. 

Since  the  flow  always  seoarates  trom  tne  centerpody 
aft  face  ana  forms  a  f ree-pressure  bounaary,  the 
reference  ideal  area  could  be  defined  as  either  of 
the  preceding  areas  minus  the  area  occupied  by  the 
separation  region,  for  example  A in  Figure  A-l. 
This  choice  has  some  merit  from  a  fluid  dynamic 
point  of  view  for  annular  nozzles,  but  it  is 
impractical  since  tne  size  of  the  separation  region 
is  difficult  to  determine.  In  the  present  study, 
the  reference  ideal  exit  area  of  the  nozzle  is 
defined  as  Aej  in  Figure  A- 1. 

Once  the  reference  ideal  choking  area  ard  the 
reference  ideal  exit  area  have  been  specified,  the 
reference  'deal  thrust  for  internal  efficiency  ano 
area  ratio  efficiency  can  be  computed  using  conven¬ 
tional  steady  one-dimensional  isentropic  nozzle 
flow  theory.  The  reference  ideal  thrust  for  pres¬ 
sure  ratio  efficiency  can  be  computed  for  a  speci¬ 
fied  value  of  nozzle  pressure  ratio  by  the  same 
techniques . 


Cg  -  A*/Ay  =  mact/(nref  ( A— 2 ) 

where  mact  is  the  actual  mass  flow  rate  and  mref 
is  the  mass  flow  rate  corresponding  to  steady 
one-dimensional  isentropic  choked  flow  at  the  geo¬ 
metric  throat  area. 


1 1 -4838 
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Fig.  A-l  Reference  areas  for  an  annular  nozzle. 


Figure  A -2  presents  the  three  nozzle  effi¬ 
ciencies  as  a  function  of  nozzle  pressure  ratio  tor 
a  nozzle  having  a  reference  ideal  area  ratio  of  3:1. 
(This  corresponds  to  an  adapted  pressure  ratio  of 
approximately  20 : 1  for  a  gas  with  a  specific  heat 
ratio  of  1.30.)  The  assumed  internal  efficiency  is 
0.97  and  the  assumed  discharge  coefficient  is  0.9A. 
The  area  ratio  efficiency  and  the  pressure  ratio 
efficiency  have  been  computed  from  conventional 
steady  one-dimensional  nozzle  flow  theory  and  the 
assumed  values  of  internal  efficiency  and  discharge 
coefficient.  The  area  ratio  efficiency  shows  the 
loss  in  performance  of  the  particular  nozzie  due  to 
the  atmospheric  pressure,  while  the  pressure  ratio 
efficiency  shows  the  loss  in  performance  due  to 
improper  matching  of  the  design  area  ratio  with  the 
operating  pressure  ratio.  The  two  efficiencies  are 
equal  only  at  the  design  area  ratio. 

The  pressure  ratio  efficiency  is  the  most 
meaningful  for  comparing  different  nozzie  desiqns 
operating  over  a  range  of  nozzle  performance  «-atios. 
Consequently,  that  efficiency  was  employed  in  the 
present  investigation. 
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Fig.  A-2  Comparison  of  nozzle  efficiency 
definitions. 
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Abstract 


An  analytical  perfomance  prediction 
methodology  for  annular  propulsive  nozzles 
with  swirl  introduced  in  the  combustor 
upstream  of  the  nozzle  is  presented.  The 
application  of  that  methodology  to  a 
specific  nozzle  design  for  a  free  vortex 
swirl  distribution  is  discussed. 
Discharge  coefficients,  specific  impulses, 
and  wall  pressure  distributions  are 
presented.  These  numerical  studies  show 
that  the  discharge  coefficient  and  the 
thrust  decrease  as  the  amount  of  swirl 
increases,  but  that  the  specific  impulse 
is  essentially  uninfluenced  by  the  swirl. 
This  methodology  will  enable  nozzle 
designers  to  account  for  the  effects  of 
swirl  in  nozzle  design. 


Nomenclature 


a  speed  of  sound 

C.-C.  constants  specifying  tangential 
1  4  velocity  distributions 

p  thrust 

m  mass  flow  rate 

M  Mach  number 

p  static  pressure 

p  stagnation  pressure 

R°  gas  constant 

S  swirl  number 

t  time 

T  static  temperature 

T  stagnation  temperature 

u°  axial  velocity  component 

v  radial  velocity  component 

w  tangential  velocity  component 

x,y  axial  and  radial  coordinate, 

respectively 

Y.  inner  radius  of  the  annular  flow 

1  geometry 

Y  outer  radius  of  the  annular  flow 

2  geometry 

y  specific  heat  ratio 

p  density 

9  tangential  coordinate 


CraHuate  Student  in  Mechanical 
Engineering.  Presently  with  Lawrence 
Livermore  Laboratory. 

+  Professor  of  Mechanical  Engineering. 
Associate  Fellow  AIAA. 

++  Professor  of  Mechanical  Engineering. 

Member  AIAA. 

Caen***  ©  »«n<«  ludtau  of  Acraoiotm  lod 
1  an  lac..  1985.  All  rttfcti  rtwrxS. 


Subscr ipts 

A  to  F  denotes  geometric  stations 

B , C , P , I VL  denote  base,  cowl,  plug, 

and  initial-value  line 
t.x.y  denotes  partial  differentiation 


Introduction 

Recent  studies  indicate  that  the 
introduction  of  swirl  ahead  of  the  combus¬ 
tor  in  axisymmetnc  dump  combustors  can 
have  very  beneficial  effects  on  the 
combustion  process.  Buckley,  Craig,  Davis 
and  Schwartzkopf 1  found  that  swirl  both 
reduced  the  reattachment  length  of  the 
combustor  flowfield  (thereby  reducing  the 
overall  combustor  length  needed  for  good 
performance),  and  helped  eliminate  des¬ 
tructive  very  low  frequency  instabilities. 
They  further  concluded  that  (in  the  range 
of  swirl  intensities  of  their  study) 
•losses  in  thrust  due  to  residual  swirl, 
at  least  to  the  sonic  point  of  the  nozzle, 
are  negligible. • 

Scharrer  and  Lilley2  made  five-hole 
pitot  probe  measurements  of  the  effects  of 
swirl  in  simulated  dump  combustors  fol¬ 
lowed  by  a  nozzle.  They  observed  a  signi¬ 
ficant  interaction  between  the  swirling 
flowfield  in  the  simulated  combustor  and 
the  nozzle  flowfield.  Their  major  objec¬ 
tive,  however,  was  the  measurement  of  the 
confined  turbulent  flow  in  the  simulated 
combustor,  not  the  nozzle  flowfield.  Con¬ 
sequently,  their  nozzles  were  simply  used 
as  downstream  blockage  components.  Both 
nozzles  used  in  their  studies  were  conven¬ 
tional  converging  nozzles  without  center- 
bodies  . 

Conley,  Hoffman  and  Thompson3 
presented  an  analytical  and  experimental 
investigation  of  the  performance  of  annu¬ 
lar  propulsive  nozzles  without  swirl.  The 
present  work  is  an  extension  of  the  per¬ 
formance  prediction  methodology  developed 
by  Conley,  Hoffman,  and  Thompson  to 
include  the  effects  of  swirl  introduced  in 
the  combustor  on  the  performance  of  annu¬ 
lar  propulsive  nozzles. 

A  similar  study  performed  by  Dutton4 
shows  trends  similar  to  the  results 
obtained  in  the  present  investigation  for 
conventional  convergent-divergent  nozzles 
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without  centerbodies  .  The  present  inves¬ 
tigation  is  concerned  with  convergent- 
divergent  nozzles  with  centerbodies . 

The  objective  of  the  present  work  was 
to  analytically  investigate  the  effect  of 
swirl  on  the  transonic  and  supersonic 
flowfields  m  annular  propulsive  nozzles, 
and  to  determine  the  effects  of  swirl  on 
mass  flow  rate,  thrust,  and  specific 
impulse . 

Performance  Prediction  Methodology- 
Geometric  Model 


The  geometric  model  considered  in 
this  investigation  is  illustrated  in  Fig. 
1.  Air  enters  at  Station  A  and  flows 
through  a  swirler  where  tangential  momen¬ 
tum  is  transferred  to  the  air  to  give  the 
desired  tangential  velocity  distribution 
at  Station  B.  Station  B  is  followed  by  a 
sudden  expansion  dump  into  the  comoustor 
inlet  at  Station  C.  Combustion  takes 
place  between  Stations  C  and  D,  where  the 
stagnation  temperature  T  rises 
corresponding  to  the  amount  of°fuel  added 
and  the  stagnation  pressure  P  descreases 
slightly  due  to  friction,  mixing,  and  heat 
addition.  The  combustion  products 
accelerate  in  the  nozzle  to  the  choked 
condition  at  the  nozzle  throat.  Station  Z, 
after  which  the  flow  continues  to 
accelerate  supersonically  from  Station  E 
to  Station  F. 


Fig.  1.  Geometric  model. 


The  concern  in  the  present  study  was 
the  effect  of  swirl  introduced  in  the 
swirler  on  the  performance  of  the  nozzle 
(i.e.,  mass  flow  rate,  thrust,  and 
specific  impulse).  Consequently,  the 
swirler,  the  sudden  expansion  dump,  and 
the  combustor  are  not  modeled  in  detail  in 
the  present  analysis.  Emphasis  is  placed 
on  the  nozzle  flowfield.  Even  so,  some 
attention  must  be  given  to  the  swirler, 
since  it  is  the  source  of  the  swirl  in  the 
nozzle  flowfield. 


air  as  it  passes  througn  the  swirler.  The 
design  of  the  swirler  and  the  interaction 
of  the  swirl  with  the  combuetion  procese 
are  of  paramount  importance  m  an  actual 
system  design.  In  the  present  etudy  of 
the  effects  of  swirl  on  nozzle  perfor¬ 
mance,  a  tangential  velocity  distribution 
representative  of  what  might  be  produced 
by  a  swirler  and  a  combustor  is  simply 
assumed  to  exist  at  the  nozzle  inlet. 

Buckely^  Craig,  Davis,  and 
Schwartzkopf  experimentally  investigated 
four  types  of  tangential  velocity  distri¬ 
butions  at  the  swirler  exit.  They  are: 


Constant  angle: 

w-C^ 

(1) 

Forced  vortex: 

w-C^y 

(2) 

Free  vortex: 

w-C3/y 

(3) 

Rankme  vortex: 

w-(C4/y) {l-exp( 

-yz/*22))U) 

where  w  is  the  tangential  velocity,  y  is 
the  radial  location,  Y.  is  the  outer 
radius  of  the  swirler,  and  C.  to  C.  are 
constants  specified  to  achieve  the  delired 
tangential  velocity  distribution.  tn 
order  to  provide  a  basis  for  comparison  of 
different  swirlers,  a  swirl  number  S  is 
employed,  which  is  defined  as 


S 


J 


2  2 
x  uwy  dy 

J  ^  u2ydy 


(5) 


The  only  tangential  velocity  distri¬ 
bution  considered  in  the  present  study  was 
the  free  vortex  distribution  given  by  Eq. 
(3).  A  centerbody  exists  throughout  the 
flowfield,  so  the  singularity  associated 
with  the  core  of  a  free  vortex  does  not 
occur  in  this  flowfield.  For  a  free  vor¬ 
tex,  yw  -  c,,  so  that  at  a  uniform  flow 
section  where  the  radial  velocity  v  »  0 
and  the  axial  velocity  u  is  uniform,  Eq. 
(5)  simplifies  to 


S 


(6) 


- - >  -  -- 

pie  function  of  the  geometry  (i.e.,  Y  ) 
and  the  axial  velocity  u.  For  a  givSn 
value  of  c3  at  the  nozzle  inlet.  Station 
Thus  1  number  SD  c»n  be  calculated. 


Swirler 


The  swirler  ie  a  set  of  axial  flow 
guide  vanes  that  Imparts  the  desired 
tangential  velocity  distribution  to  the 


SD  ' 


D  2D 


(7) 
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The  corresponding  values  of  swirl 
number  at  the  combustor  inlet  S  ,  swirler 
exit  S„,  and  nozzle  throat  S_,  can  be 
estimated  from  the  known  value  of  Sn  and 
Eq .  (6).  The  major  process  occuring 
between  Stations  C  and  D  is  the  increase 
in  stagnation  temperature  T  due  to 
combustion.  For  a  Rayleigh0  line  flow 
(i.e.,  heat  addition  at  constant  area)  m 
the  absence  of  swirl,  it  can  be  shown  that 
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Thus,  for  known  values  of  u  ,  T  , 
ToC'  UC  can  fae  calculated,  so°that° 


(8) 


and 


SD  _  ^3/UDY2D  _ 

C3/“ci:2D  UD 


(9) 


by  the  nozzle  la  directly  related  to  the 
nozzle  exhaust  velocity.  For  rocket 
engines,  the  gross  thrust  produced  by  the 
nozzle  is  al3o  the  net  thrust.  However, 
for  airDreathing  engines  (turbines  and 
ramjets),  the  net  thrust  is  the  difference 
between  the  nozzle  grosB  thrust  and  the 
inlet  ram  drag.  Consequently,  all  losses 
in  nozzle  gross  thrust  for  airbreathers 
are  magnified  several  times  m  net  thrust, 
and  it  is  imperative  that  the  nozzle  effi¬ 
ciency  be  as  large  as  possible. 

Figure  2  illustrates  the  thrust  pro¬ 
ducing  components  of  a  generic  annular 
nozzle.  The  nozzle  thrust  is  composed  of 
the  thrust  developed  across  the  supersonic 
initial-value  line  F  ,  the  thrust 
developd  by  the  pressure  acting  on  the 
cowl  P  ,  the  thrust  developed  by  the  pres¬ 
sure  acting  on  the  plug  F  ,  and  the  thrust 

acting  on  tne  plug  base  F  . 
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At  the  sudden  expansion  dump  loca¬ 
tion,  Yj  changes  from.  Y^-  to  *2C'  FoE 
constant  mass  flow  rate  m,  and  neglecting 
compressibility  effects,  which  are  small, 
u  is  inversely  proportional  to  the  flow 
area,  which  is  directly  proportional  to 
Y^.  Thus, 


2C'  2B 1 


(10) 


Substituting  this  result  into  Eq.  (6) 
yields 


Fig.  2.  Thrust  producing  components  of 
an  annular  nozzle. 
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At  the  nozzle  throat,  Station  E,  u 
is  three  to  five  times  as  large  as  un,  anS 
Y2E  may  be  an>rwh® re  from  50  to  90  percent 
of  Y^  .  When  specific  values  are  given, 

Sp  can  be  calculated  from  S„. 

D 

Consequently,  for  a  specified  value 
of  the  swirl  yw  at  the  nozzle  inlet,  the 
swirl  number  SQ  at  the  nozzle  inlet  can  be 
calculated.  Knowing  the  other  motor 
operating  parameters,  tne  values  of  the 
swjri  number  S  can  be  estimated  at  other 
locations  m  the  motor. 


The  relative  magnitudes  of  these 
thrust  components  depend  on  the  nozzle 
design.  For  nozzie  geometric  area  ratios 
between  3  and  5,  as  considered  m  this 
study,  the  magnitude  of  the  supersonic 
initial- value  line  thrust  is  approximately 
80  to  85  percent  of  the  total  thrust,  the 
cowl  thrust  is  from  10  to  15  percent  of 
the  total  thrust,  and  the  plug  and  base 
thrusts  comprise  the  remainder. 

The  supersonic  initial-value  line 
thrust  F[Vt  is  given  by 

P  i  vl  ”  ;vl  u  dm  *  ivl  P  27ry  dy 


Performance  Model 


The  cowl  thrust  F^  is  given  by 


The  propulsive  nozzle  is  the  thrust 
producing  component  of  jet  propulsion 
engines  such  as  rockets,  turbojets  and 
ramjets.  The  function  of  the  propulsive 
nozzle  is  to  convert  the  random  thermal 
energy  of  the  gases  entering  the  nozzle 
into  directed  kinetic  energy  at  the  nozzle 
exit,  thus  creating  a  thrust  on  the  pro¬ 
pulsive  system.  The  expansion  process 
must  be  accomplished  as  efficiently  as 
possible,  since  the  gross  thrust  developed 


FC  *  Cowl  P  dy 

The  plug  thrust  Fp  is  given  by 

PP  “  Plug  P  2,7y  dy 


(13) 


(14) 
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The  plug  base  thrust  Fg  is  given  by 


where  P  is  the  average  base  pressure  act¬ 
ing  on  the  truncated  plug  case  area  Ag . 

From  Eas.  (12)  to  (15),  it  is  seen 
that  the  flow  properties  across  the  super¬ 
sonic  init lal-value  line,  the  pressure 
distributions  along  the  cowl  and  the  plug, 
and  the  base  pressure  must  be  calculated 
to  enable  the  calculation  of  tne  nozzle 
performance. 

Flowfield  Model 


The  VNAP  code  was  U3ed  as  the 
auosomc/transonic  flowfield  prediction 
code  for  the  annular  nozzle  performance 
prediction  methodology  developed  by  Con¬ 
ley,  Hoffman,  and  Thompson3.  The  VNAP 
code  solves  the  Nav let -Stokes  equations 
for  two-dimensional  (planar  or  axisym- 
metric),  time  dependent,  compressible  flow 
usina  the  second-order  accurate,  MacCor- 
mack  finite  difference  scheme.  lnviscid 
flowfields  can  be  >ived  by  setting  all  of 
the  viscous  teri.*'  to  zero,  wmch  yields 
the  Euler  equations.  An  explicit  artifi¬ 
cial  viscosity  is  included  for  shock  wave 
calculations.  The  fluid  is  assumed  to  be 
a  perfect  gas.  The  steady  state  solution 
is  obtained  as  the  asymptotic  solution  for 
large  time. 


The  accuracy  and  efficiency  of  nozzle 
performance  prediction  techniques  are 
highly  dependent  on  the  flowfield  model 
chosen  to  represent  the  actual  flowfield. 
In  the  present  study,  the  flowfield  modei 
is  based  on  the  following  assumptions: 

1.  steady  axisymmetric  flow, 

2.  mviscid  nonconducting  fluid, 

3.  no  body  forces, 

4.  thermally  and  caloricaily  perfect  gas, 
and 

5.  separation  and  base  pressures  pre¬ 
dicted  from  empirical  correlations. 

The  governing  equations  consist  of  the 
continuity  equation,  the  component  momen¬ 
tum  equations,  the  energy  equation  and  the 
thermal  and  caloric  equations  of  state.  A 
detailed  discussion  of  these  equations  is 
presented  in  Ref.  (5). 

The  choice  of  a  numerical  method.  or 
methods,  to  solve  these  governing  equa¬ 
tions  is  also  crucial  to  the  success  of 
the  performance  prediction  procedure.  A 
sub3onic/transonic  solution  is  required  to 
define  the  supersonic  initial-value  line 
thrust,  while  a  supersonic  solution  is 
required  to  define  the.  cowl,  plug,  and 
base  thrusts.  Experience J  has  shown  that 
subsonic/tranaonic  flowfields  can  be  ade¬ 
quately  determined  from  the  steady  state 
solution  to  finite  difference  time  march¬ 
ing  solutions,  whereas  the  method  of 
characteristics  is  the  most  accurate  means 
for  determining  supersonic  flowfields  in 
nozzles.  Therefore,  the  numerical  methods 
selected  for  the  performance  prediction 
procedure  in  the  present  study  are  a  fin¬ 
ite  difference  time  marching  approach  to 
define  the  subsonic/transonic  flowfield, 
and  a  method  of  characteristics  approach 
to  define  the  supersonic  flowfield. 

Subsonic/Transonic  Flowfield 

The  finite  difference  code  developed 
to  solve  the  aubsomc/tr ansonic  flowfield 
is  a  modification  of  the  VNAP  (Viscous 
Nozzle  Analysis  Program)  code  developed  by 
Cline6  to  include  the  effects  of  swirl. 
The  resulting  code  is  called  SNAP  (Swir¬ 
ling  Nozzle  Analysis  Program) . 


The  computat lonal  grid  employed  by 
the  VNAP  code  is  a  transformed  equally 
spaced  grid.  Equally  spaced  lines  of  con¬ 
stant  £  normal  to  the  x-axis  comprise  one 
of  the  transformed  coordinates.  Each 
vertical  line  segment  between  the  plug  and 
the  cowl  is  divided  into  a  number  of 
equally  spaced  increments.  Lines  of  con¬ 
stant  7f  join  these  points.  The  physical 
arid  is  illustrated  in  Fig. 
transformation  eauations  are 


C  -  x 


and 


y  -  yc(x) 
yc(x;  -  yp;x) 


is  the  plug  contour. 


3  . 

The 

(16) 

TT 

and 

ypU) 

Fig.  3-  Typical  SNAP  grid. 


In  the  present  investigation,  the 
VNAP  code  was  modified  to  account  for  the 
presence  of  3wirl  m  the  flowfield  by 
add inq  the  tangential  momentum  equation  to 
the  set  of  governing  equations  and  the 
tangential  velocity  component  to  all  the 
aovernma  eauations.  where  appropriate. 
Since  the  flowfield  is  axisymmetric,  all 
derivatives  with  respect  to  the  tangential 
direction  0  are  zero,  even  though  the 
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tangential  velocity  component  w  is 
nonzero.  The  governing  equations  i or  this 
flowfield  are: 


pt  +  upx  +  Vpy  +  pUx  +  pvy  +  pV/y  "  0 
put  +  puux  +  pvuv  +  Px  -  0 

pvt  +  puvx  +  pvvy  +  Py  -  pw2/y  -  0 


(1?) 


(18  ) 


(19) 


+  uw  + 
t  X 


vw  + 

y 


vw/y 


0 


(20) 


P  +  uP  +  vP  -  a2(p  +up  +vp  i  -  0  (21) 
u  x  y  u  x  y 

where  u,  v,  and  w  are  the  velocity  com¬ 
ponents  in  the  x,  v,  and  6  directions, 
respectively,  P  is  the  pressure,  p  is  the 
density,  a  is  the  speed  of  sound,  and  sub¬ 
scripts  denote  partial  differentiation. 
For  a  perfect  gas. 


P  -  pRT 


(22) 


a2  -  yP/p 


(23) 


Equation  (20)  can  be  written  as 


D(yw) 

Dt 


0 


(24) 


which  is  in  substantial  derivative  form. 
Equation  (24)  states  that  the  swirl  yw 
remains  constant  on  the  pathlines  of  an 
unsteady  flowfield  and  on  the  streamlines 
of  a  steady  flowfield.  Consequently,  if 
the  distribution  of  the  swirl  yw  is  speci¬ 
fied  at  the  nozzle  inlet,  that  swirl  dis¬ 
tribution  will  be  preserved  throughout  the 
flowfield.  Thus,  two  modifications  to 
VNAP  were  made  to  include  the  effects  oi 
swirl:  (1)  tracking  the  streamlines  of 
the  flowfield  so  that  the  values  of  swirl 
yw  at  the  nozzle  inlet  can  be  tracked 
through  the  flowfield  as  the  streamlines 
change  radial  position,  and  (2)  adding  the 
term  pw  /y  to  Eq.  (19),  where  w  is  com¬ 
puted  from  the  known  value  of  swirl  yw 
along  each  streamline. 

The  initial  tangential  velocity  dis¬ 
tribution  w(y)  is  specified  at  Station  B, 
illustrated  in  Fig.  4.  Station  B  is  an 
upstream  station  ahead  of  the  nozzle 
inlet.  Typically,  Station  B  is  the  exit 
of  tbe  swirler.  The  mass  flow  distribu¬ 
tion  m(y)  is  computed  across  Station  B,  so 
that  the  swirl  yw  is  a  known  function  of 
A.  This  function,  yw(m),  is  then  U3ed  to 
compute  the  swirl  yw  at  each  point  in  the 


Fig.  4.  Swirl  specification  model. 


nozzle  m  the  following  manner.  After 
each  time  step,  the  mass  flow  rate  is 
integrated  across  the  flowfield  at  each 
axiai  location,  from  the  centerboay 
(m  -  0)  to  the  nozzle  wail  (m  -  "'total ^  ' 
The  fractional  mass  flow  rate  at  each 
radial  station  is  calculated  by  dividing 
the  calculated  mass  flow  rate  at  that 
radial  station  by  the  total  mass  flow  rate 
at  that  axial  location.  The  resulting 
fractional  mass  flow  rate  is  then  U3ed  to 
interpolate  for  the  swirl  yw  at  the 
corresponding  computational  grid  point  by 
quadratic  interpolation  of  the  swirl  func¬ 
tion  yw(m)  computed  at  Station  B.  The 
swirl  velocity  w  is  then  computed  by 
dividing  the  interpolated  value  of  the 
swirl  yw  by  the  value  of  y  at  the  computa¬ 
tional  grid  point. 

This  procedure  is  consistent  with  the 
governing  partial  differential  equations 
only  at  the  steady  flow  limit,  because 
actual  unsteady  pathlines  are  not  tracked. 
Essentially,  this  procedure  matches  a  nor¬ 
malized  steady  mass  flow  rate  at  the 
swirler  exit  to  a  normalized  mass  flow 
rate  at  each  point  in  the  nozzle.  In  an 
unsteady  flow,  these  two  points  would  not 
necessarily  lie  on  the  same  pathline. 
However,  at  the  steady  flow  limit,  these 
two  points  do  lie  on  the  same  streamline. 

In  a  typical  run,  800  to  1000  time 
steps  are  taken  by  the  SNAP  code  to  relax 
the  initially  assumed  one-dimensional 
flowfield  to  a  steady  state.  An  initial- 
value  line  along  which  the  projection  of 
the  Mach  number  in  the  xy  plane  is  super¬ 
sonic  is  then  specified  across  the  throat 
region.  Values  of  the  projected  Mach 
number  and  flow  angle  along  that  line  are 
determined  from  the  SNAP  flowfield  by 
interpolation.  Least  squares  quadratic 
bivariate  interpolating  polynomials  are 
fit  to  the  nine  closest  SNAP  grid  points 
to  each  initial-value  line  point.  In  this 
manner,  a  supersonic  initial-value  line  is 
determined,  which  can  be  used  to  initiate 
the  solution  of  the  supersonic  flowfield 
by  the  method  of  characteristics. 


5 


Supersonic  Flowfield 


The  supersonic  flowfield  is  calcu¬ 
lated  by  the  method  of  character istics  for 
steady  two-dimensional  flow.  For  steady 
flow,  Eqs.  (17)  to  (21)  become 


upx  +  vpy  +  pu^  +  pvy  +  pv/y  -  0  (25) 


puux  +  pvuy  +  Px  -  0  (26) 

puv^  +  pWy  +  -  pw2/y  -  0  (27) 

uw  +■  vw  +  vw/y  -  0  (28) 

*  y 

uPx  +  vPy  "  a2(upx  +  vPy)  “  o  (29) 


Equation  (28)  can  be  written  as 

^  •  0 


Fig.  5.  Typical  left-running  Mach  lines 
for  supersonic  annular  swirling 
f  low . 

continued  until  the  nozzle  exit  lip  point, 
point  E.  is  determined.  Left-running  Mach 
lines  are  then  terminated  at  the  last 
right-running  Mach  line  emanating  from 
Point  E.  This  procedure  is  continued 
until  point  F,  at  the  end  of  the  plug  is 
determined.  The  base  pressure  on  the  face 
of  the  truncated  plug  is  computed  from  the 
following  empirical  correlation: 


which  states  that  the  swirl  yw  is  constant 
along  streamlines. 


base 


0.846  PF/M^'3 


<31; 


The  derivation  of  the  characteristic 
and  compatibility  equations  corresponding 
to  Eqs.  (25)  to  (29)  is  a  straightforward, 
although  somewhat  lengthy,  procedure.  The 
derivation  is  presented  by  Kornblum  and 
Thompson  The  results  are  not  presented 
here.  In  brief,  three  characteristics  are 
determined  m  the  xy  plane:  the  stream¬ 
line  and  the  right-  and  left-running  Mach 
lines.  Three  compatibility  equations  are 
determined  along  the  streamline  and  one 
compatibility  equation  is  determined  along 
each  Mach  line.  These  five  compatibility 
equations  are  used  to  solve  for  u,  v,  w, 
P ,  and  p . 

A  direct  marching  method  of  charac¬ 
teristics  procedure  is  employed  to  solve 
the  characteristic  and  compatibilty  equa¬ 
tions.  The  solution  procedure  constructs 
the  Mach  line  network  illustrated  in  Fig. 
S.  The  characteristic  and  compatibility 
equations  are  integrated  numerically  by 
the  second-order  accurate  modif ied-Euler 
predictor-corrector  method.  The  solution 
is  initiated  from  the  supersonic  initial- 
value  line,  line  TT‘ ,  determined  from  the 
SNAP  analysis.  Left-running  Mach  lines 
are  emanated  from  line  TT’,  starting  at 
point  T,  and  continued  across  the  flow- 
field  until  they  intersect  the  cowl,  con¬ 
tour  TE.  The  last  such  Mach  line  from 
point  T*  intersects  the  cowl  at  point  T" . 
Left-running  Mach  lines  are  then  ori¬ 
ginated  from  points  on  the  centerbody  and 
continued  across  the  flowfield  to  inter¬ 
sect  the  nozzle  wall.  This  procedure  is 


The  above  procedure  determines  all  of 
the  data  required  to  evaluate  the  thrust 
components  illustrated  in  Fig.  2. 

Resu  Its 

Nozzle  flowfields  were  computed  for 
12  cases.  Figure  6  defines  the  flowfield 
geometry  for  all  12  cases.  Only  free  vor¬ 
tex  swirlers  were  considered.  The  swirl 
strength  was  characterized  by  the  yw  pro¬ 
duct,  where  y  is  the  flow  radius  (ft),  and 
w  is  the  trangenital  velocity  component 
(ft/sec).  Values  of  the  swirl  yw  of  0, 
50,  100,  150,  200,  and  250  ftZ/sec  were 
analyzed  for  stagnation  temperatures  at 
the  nozzle  inlet  of  2500  R  and  3500  R. 
The  stagnation  pressure  was  35  psia  in 
every  case.  In  all  cases,  the  specific 
heat  ratio  y  -  1.40  and  the  gas  constant  R 
-  53.35  (f t-lbf ) / ( lbm-R) .  The  scope  of 
the  parametric  studv  is  defined  in  Table 
1. 
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Table  1.  Scope  of  the  parametric  stud/. 


Run 

Swirl  yw  T 

0 

(ft^/secj  R 

1 

0 

2500 

2 

50 

2500 

3 

100 

2500 

4 

150 

2500 

5 

200 

2500 

6 

250 

2S00 

7 

0 

3500 

8 

50 

3500 

9 

100 

3500 

10 

150 

3500 

11 

200 

3  500 

12 

250 

3500 

The  swirl  number  S  defined  by  Eq.  (5) 
varies  with  axial  location  in  the  flow- 
field.  At  the  nozzle  inlet  (x  -  -  6.0  in. 
foi;  the  SNAP  analysis),  a  swirl  yw  of  50 
ft  /sec  corresponds  approximately  to  a 
swirl  number  S  -  0.14  for  T  -  2500  R,  and 
S  -  0.12  for  Tq  -  3500  8.  The  swirl 

number  S  varies  approximately  linearly 
with  the  swirl  yw  at  that  location,  so 
that  the  range  of  yw  values  considered  in 
this  study  (i.e.,  0  to  250  ft^/sec) 

corresponds  to  swirl  numbers  between  0  and 
0.7  at  the  nozzle  inle*-  At  the  nozzle 
throat,  the  swirl  numbers  are,  of  course, 
much  smaller.  For  the  nozzle  geometry 
considered  here,  the  swirl  numbers  at  the 
throat  are  about  30  percent  of  the  values 
at  the  nozzle  inlet. 

The  swirl  number  upstream  of  the 
combustor  depends  on  the  amount  of  heat 
added  in  the  combustor,  as  indicated  by 
Eqs.  (8)  and  (9).  Por  stagnation  tempera¬ 
tures  of  approximately  1000  R  entering  the 
combustor  and  a  swirl  yw  -  50  ft  /sec,  S  « 
0.37  for  T  -  2500  R  and  S  -  0.44  for  T 
3500  R.  The  swirl  number  varies  approxi¬ 
mately  linearly  with  yw,  so  the  range  of 
yw  values  considered  in  this  "study 
corresponds  to  swirl  numbers  between  0  and 
2.20  at  the  combustor  inlet. 

The  swirl  number  at  the  swirler  exit, 
upstream  of  the  sudden  expansion  dump, 
depends  on  the  geometry  of  the  sudden 
expansion,-  as  expressed  by  Eq.  (11).  For 
yw  -  50  f t vsec ,  S  -  0.24  for  T  -  2500  R 
at  the  nozzle  inlet,  and  S  -  0.?9  for  T 
3S00  R  at  the  nozzle  inlet.  The  sw?rl 
number  varies  approximately  linearly  with 
yw,  so  the  range  of  yw  values  considered 
in  this  study  corresponds  to  swirl  numbers 
between  0  and  1.45  at  the  swirler  exit. 

Figure  7  presents  constant  Mach 
number  lines  in  the  transonic  region  for 
Ca3e  7  (T  -  3S00  R,  yw  »  0)  as  computed 
by  the  SNAP  code.  Comparisons  with  the 
constant  Mach  number  lines  for  the  other 
11  cases  show  differences  in  the  third 
significant  figure  for  most  quantities. 


Pig.  7.  Transonic  Mach  number  distribution 
f or  Case  7 . 

Thus,  Fig.  7  is  representative  of  the 
plots  for  all  12  cases.  There  are  some 
differences,  however,  which  are  discerni¬ 
ble  m  the  integrated  mass  flow  rate 
values  tabulated  in  Table  2. 

Figures  8  and  9  present  the  left- 
running  Mach  line  pattern  in  the  xy-plane 
for  the  supersonic  flowfields  for  Cases  1 
and  6,  respectively.  Case  1  is  the  no 
swirl  c^se  with  T  -  2500  R.  Case  6 ,  yw  - 
250  ft  /sec  and°  T  -  2500  R,  represents 
the  largest  swirl  case  considered  in  the 
study.  In  Figs.  8  and  9,  the  left-running 
Mach  lines  emanating  from  the  plug  surface 
are  terminated  at  the  right-running  Mach 
line  emanating  from  the  end  of  the  cowl. 
Comparison  of  Figs.  8  and  9  shows  that 
swirl  has  the  effect  of  compressing  the 
flow  on  the  plug  surface,  resulting  m 
coalescing  left-running  Mach  lines. 


i 


x,  in. 


Fig*  8.  Left-running  Mach  lines  for 
Case  1. 


Fig.  9  -  Left-running  Mach  lines  for 
Case  6. 


Figure  10  is  a  comparison  of  the 
pressure  distributions  on  the  plug  for 

Cases  1  and  6.  There  is  a  substantial 
pressure  increase  on  the  plug  due  to  the 
effects  of  swirl.  Unfortunately,  the  area 
over  which  the  pressure  increase  acts  is 
small,  so  the  increase  m  thrust  is  only 
modest.  The  pressure  distributions  for 
the  intermediate  swirl  intensities  He 
between  the  results  of  Cases  1  and  6  in  a 
very  orderly  fashion. 


0  1  2  3  4  5 

«,  la 


Fig.  10.  Pressure  distributions  on  the 
plug  for  Cases  1  and  6. 


Pigure  11  presents  the  pressure  dis¬ 
tribution  on  the  cowl  for  Cases  1  and  6. 
The  pressure  distributions  on  the  cowl 
surface  are  almost  indistinguishable  for 
all  swirl  values.  There  are  some  differ¬ 
ences,  nowever,  which  are  discernable  m 
the  integrated  thrust  values  tabulated  in 
Table  2.  Because  the  cowl  projected  area 
is  much  larger  than  the  plug  area,  a  much 
smaller  pressure  difference  is  required  to 
produce  distinguishable  differences  in  the 
integrated  thrust  values  on  the  cowl. 


Fig.  11.  Pressure  distributions  on  the 
cowl  for  Cases  1  and  6. 


Table  2  summarizes  the  mass  flow 
rates,  thrusts,  and  specific  impulses  for 
the  12  cases  studied.  Cases  1  to  6  are 
for  To  -  2500  R,  and  Cases  7  to  12  are  for 
T  -  3500  R.  The  reference  (one- 

dlmensional)  mass  flow  rate  for  Cases  1  to 
&  iS  ml-D  *  21.9842  lbm/sec,  and  for  Cases 
J  to  lZ,  m,  “  18.5800  lbm/sec.  Mass 
flow  values  were  determined  by  integrating 
across  the  initial-value  line  that  is  gen¬ 
erated  from  the  results  of  the  SNAP 
analysis.  The  discharge  coefficients  are 
plotted  in  Fig.  12.  The  values  decrease 
with  increasing  swirl  in  a  smooth,  con¬ 
sistent  manner  and  indicate  consistency  in 
the  initial-value  data  obtained  from  the 
SNAP  analysis.  Thrust  values  were  com¬ 
puted  from  the  results  of  the  method  of 
characteristics  (MOC)  analysis.  The 
individual  thrust  contributions  from  the 
initial-value  line,  the  cowl,  the  plug, 
and  the  base  (illustrated  In  Fig.  2)  are 
tabulated  separately.  All  values  are 
referenced  to  zero  ambient  pressure.  The 
base  pressures  were  computed  from  Eq 
(31). 


Table  2.  Comparison  of  performance  parameters. 


Run 

Swirl  yw  Tq 

(ft2/sec)  R 

Mass  Flow 

( lbm/sec ) 

CD 

I 

sp 

{  lbf -sec/lbm) 

1 

0 

2500 

21.491 

0.9776 

141. 594 

2 

50 

2  son 

21.456 

0.9760 

141.625 

3 

100 

2500 

21.350 

0.9712 

141.665 

4 

150 

2500 

21.173 

0.9631 

141.753 

5 

200 

2500 

20.929 

0.9520 

141.894 

6 

250 

2500 

20.611 

0.9375 

142.161 

7 

0 

3500 

18.163 

0.9776 

167.538 

8 

50 

3500 

18.142 

0.9764 

167.569 

9 

100 

3500 

18.078 

0.9730 

167.602 

10 

150 

3500 

17.971 

0.9672 

167.670 

11 

200 

3500 

17.823 

0.9593 

167.774 

12 

250 

3500 

17.633 

0.9490 

167.952 

Run 

1  v  1 

Cowl 

Thrust  (lbf) 

Plug 

Base 

Total 

1 

2582.08 

298.26 

155.34 

7.32 

3043.00 

2 

2577.11 

298.06 

155.57 

7 . 96 

3038.70 

3 

2562.23 

297.45 

156.28 

8.59 

3024.55 

4 

2537.58 

296.44 

157.53 

9.  79 

3001.34 

5 

2503.38 

295.02 

159.37 

11.93 

2969.70 

6 

2459.43 

292.67 

162.17 

15.82 

2930.09 

7 

2582.08 

298.26 

155.34 

7.32 

3043.00 

8 

2578.53 

298.11 

155.51 

7.90 

3040.65 

9 

2567.89 

297.68 

156.01 

8.34 

3029.92 

10 

2550.23 

296.96 

156.88 

9.14 

3013.21 

11 

2525.67 

295.95 

158.17 

10.45 

2990.24 

12 

2494.37 

294.63 

159.87 

12.63 

2961.50 

Fig.  12.  Discharge  coefficients. 


The  following  observations  can  De 
made  from  the  results  obtained  in  this 
study. 

1.  The  data  are  extremely  consistent 
giving  a  high  degree  of  confidence  m  the 
trends.  (We  also  believe  the  absolute 
values  are  reliable.)  The  thrust  values 
for  Cases  1  and  7  are  essentially  identi¬ 
cal.  This  is,  of  course,  what  one  would 
hope,  since  changing  the  stagnation  tem¬ 
perature  does  not  affect  the  absolute 
thrust  values.  However,  the  initial-value 
line  calculated  by  the  TNAP  code  might  be 
expected  to  be  somewhat  different,  sines 
the  flow  speeds  and  permissible  time  steps 
for  those  calculations  are  different. 
This  result  gives  confidence  in  the  con¬ 
sistency  of  the  results. 

2.  For  t,hose  applications  for  which 
yw  <  100  ftVsec.  the  effect  of  swirl  on 
the  mass  flow  rate  and  thrust  values  is 
very  small  and  can  probably  be  ignored  in 
preliminary  design  and  overall  performance 
considerations-  It  is  probable  that  this 
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observation  is  ntu  valid  for  oinai  than 
free  vortex  swiriers. 

3.  Even  at  hignet  awni  values.  me 
specific  impulse  of  me  nczzie  ia  not 
degraded  dv  swirl  The  calculations  snow 
some  increase  in  1  with  swirl.  There¬ 
fore,  if  high  swirl  is  to  he  used.  me 
effects  of  swirl  on  nozzle  t.nrust  can  oe 
compensated  for  ;n  tne  design  dv  increas¬ 
ing  the  throat  area  to  compensate  for  me 
decreased  mass  flow  rate,  thus  increasing 
the  thrust  to  tne  required  value  without 
any  loss  m  nozzle  efficiency. 


4 .  When 
are  contoured 
the  optimum 
affected  dv 
f lowf le Id . 


the  cowl  and  plug  surfaces 
to  optimize  tne  performance, 
contour  will  probably  oe 
the  presence  of  swirl  in  the 


Cone  1  us  l or.s 


An  analytical  performance  prediction 
methodology  for  analyzing  tne  effects  ei 
swirl  in  annular  propulsive  nozzles  nas 
Deer,  developed.  Computed  :  eauits  for  a 
aeries  of  swirl  intensities  for  tree  vor¬ 
tex  swiriers  are  extremely  consistent  and 
indicate  a  nigh  reliability  m  tne 
results.  Swir„  decreases  tne  disenarge 
coeffi-ient  of  a  nozzle  Dut  nas  very  lit¬ 
tle  effect  on  tne  vacuum  specific  impulse. 

For  swirl  intensities  often  encoun¬ 
tered  in  rarrnet  and  turnout  applications, 
the  effect  of  swirl  or.  the  nozzle  perfor- 
mace  is  small,  and  can  prooaoiy  oe 
neglected  m  preliminary  sea ign  and  per¬ 
formance  calculations. 

Experimental  verification  of  tne 
results  is  needed.  Application  of  tne 
methodology  to  other  than  free  vortex 
swiriers  is  needed. 
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Abstract 

An  analytical  performance  prediction  methodology  for 
annular  propulsive  nozzles  with  swirl  introducer,  in  the 
combustor  upstream  of  the  nozzle  is  presented  The 
methodology  is  applied  to  investigate  the  effects  of  swirler 
design  on  the  performance  of  annular  propulsive  nozzles. 
Four  types  of  swirlcrs  were  investigated  free  vortex, 
constant  angle,  forced  vortex,  and  Rankine  vortex  swirlcrs. 
Discharge  coefficients,  specific  impulses,  and  wall  pressure 
distributions  are  presented.  These  numerical  studies  show 
that  the  discharge  coefficient,  the  thrust,  and  the  vacuum 
specific  impulse  decrease  as  the  amount  of  swirl  increases, 
but  that  the  decrease  in  specific  impulse  is  modest.  The 
effect  of  swirl  on  the  discharge  coefficient,  for  all  four 
swirler  designs,  correlates  well  with  the  mass-averaged  swirl 
introduced  into  the  flow  field  by  the  sw  irler.  However,  the 
decrease  in  vacuum  specific  impulse  is  a  function  of  the 
swirler  design.  The  forced  vortex  swirler  has  the  least 
decrease  in  specific  impulse  witn  increasing  swirl  while  the 
free  vortex  swirler  has  the  greatest  decrease.  This 
methodology  will  enable  nozzle  designers  to  account  for  the 
effects  of  swirl  in  nozzle  design. 


Nomenclature 


a  speed  of  sound 

CpC*  constants  specifying  tangential  velocity 

distributions 

Cj-Cg  constants  used  to  set  the  mass  flow  rate 
F  thrust 

m  mass  flow  rate 

M  Mach  number 

P,P0  static  and  stagnation  pressure,  respectively 
R  gas  constant 

S  swirl  number 

t  time 

T,T0  static  and  stagnation  temperature,  respectively 

u,v,w  velocity  components 

V  velocity  magnitude 

x.y  axial  and  radial  coordinates,  respectively 

vw  mass-averaged  sw  ir! 

Yj.Yz  inner  and  outer  radii  of  the  flow  passage, 
respectively 

y  specific  heat  ratio 

p  density 

8  tangential  coordinate 

n,  t)i  transformed  coordinates 
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Subscripts 


A  to  F  geometric  stations 

B.C.P.1VL  base,  cowl,  plug,  and  initial-value  line, 
respectively 


Introduction 

Recent  studies  indicate  that  the  introduction  of  swirl  ahead 
of  the  combustor  in  axisymmetric  dump  combustors  can 
have  very  beneficial  effects  on  the  combusnon  process. 
Buckley,  Craig,  Davis  and  Schwartzkopf1  found  that  swirl 
both  reduced  the  reattachment  length  of  the  combustor 
flowfield  (thereby  reducing  the  overall  combustor  length 
needed  for  good  performance),  and  helped  eliminate 
destructive  very  low  frequency  instabilities.  They  further 
concluded  that  (in  the  range  of  swirl  intensities  of  their 
study)  "losses  in  thrust  due  to  residual  swirl,  at  least  to  the 
sonic  point  of  the  nozzle,  are  negligible." 

Scharrer  and  Lilley2  made  five-hole  pitot  probe 
measurements  of  the  effects  of  swirl  in  simulated  dump 
combustors  followed  by  a  nozzle.  They  observed  a 
significant  interaction  between  the  swirling  flowfield  in  the 
simulated  combustor  and  the  nozzle  flowfield.  Their  major 
objective,  however,  was  the  measurement  of  the  confined 
turbulent  flow  in  the  simulated  combustor,  not  the  nozzle 
flow  field.  Consequently,  their  nozzles  were  simply  used  as 
downstream  blockage  components.  Both  nozzles  used  in 
thi’r  studies  were  conventional  converging  nozzles  without 
centerbodies. 

Conley,  Hoffman  and  Thompson3  presented  an  analytical 
and  experimental  investigation  of  the  performance  of 
annular  propulsive  nozzles  without  swirl.  Kornblum, 
Thompson,  and  Hoffman4  extended  the  performance 
prediction  methodology  developed  by  Conley  et  ai  to 
include  the  effects  of  swirl  introduced  in  the  combustor  on 
the  performance  of  annular  propulsive  nozzles  for  free 
vortex  swirlcrs  A  similar  study  performed  by  Dutton' 
shows  trends  similar  to  the  results  obtained  by  Kornblum.  et 
al  for  conventional  converging-diverging  nozzles  without 
centerbodies 

The  present  work  is  an  extension  of  the  performance 
prediction  methodology  developed  by  Kornblum. 
Thompson,  and  Hoffman4  to  include  constant  angle 
swirlers.  forced  vortex  swirlers,  and  Kankme  vortex 
swirlers 


Performance  Prediction  Methodology 
Geometric  Model 

The  geometric  model  considered  in  this  investigation  is 
presented  in  Fig  1.  Air  enters  at  Station  A  and  flows 
through  a  swirler  where  tangential  momentum  is  transferred 
to  the  air  to  give  the  desired  tangential  velocity  distribution 
at  Station  B  Station  B  is  followed  b\  a  sudden  expansion 
dump  into  the  combustor  inlet  at  Station  C  Combustion 
takes  place  between  Stations  C  and  D.  where  the  stagnation 
temperature  T0  rises  corresponding  to  the  amount  of  fuel 
added,  and  the  stagnation  pressure  P„  decreases  slightly  due 
to  friction,  mixing,  and  heat  addition.  The  combustion 
products  accelerate  in  the  nozzle  to  the  choked  condition  at 
the  nozzle  throat.  Station  E,  after  which  the  flow  continues 
to  accelerate  supersonically  from  Station  E  to  Station  F 

The  concern  in  the  present  study  was  the  effect  of  swirl 
induced  in  the  swirler  on  the  performance  of  the  nozzle 
fi.e.,  mass  flow  rate,  thrust,  and  specific  impulse) 
Consequently,  the  swirler,  the  sudden  expansion  dump,  and 
the  combustor  are  not  modeled  in  detail  in  the  present 
analysis  Emphasis  is  placed  on  the  nozzle  flouficld  The 
swirJ  introduced  by  the  swirler  is  assumed  to  flow  through 
the  combustor  and  nozzle  unchanged  in  magnitude 


which  is  the  substantial  derivauve  form  of  the  equation 
For  steads  How.  Eq  (5)  requires  that  the  swiri  yw  remain 
constant  on  a  streamline  Thus,  the  swirl,  yw,  as  a  function 
of  the  stream  function,  is  the  same  throughout  the  flow  field 
Irom  Station  B  (the  swirler  exit)  to  Station  F  (the  nozzle 
exit;  Furthermore,  the  mass-averaged  swirl.  Fu.  is  me 
same  at  all  stations  from  the  swirler  exit  to  the  nozzle  exit, 
where 


yw 


|  1  -  j  1  « 

—  f  vw  dm  =  —  f .  vw  pu  \  ds 

m  V  "  n»  V  ' 


(b, 


The  swirl  distribution  at  the  nozzle  inlet.  Station  D.  is 
determined  as  follows  The  swirl  distribution  at  the  swirler 
exit.  Station  B.  is  calculated  as  a  function  of  the  fractional 
mass  flow  rate  That  sane  swirl  distribution  is  specified  at 
the  nozzle  entrance.  Station  D,  as  initial  conditions  for  the 
transonic  flow  analysis  The  swirl  distribution  at  Station  B 
is  computed  as  a  function  of  the  fractional  mass  flow  rate  as 
follow  s 


(a)  The  stagnation  temperature  and  pressure  at 
the  swirler  inlet.  Station  A.  are  assigned 


Performance  Model 


(b)  The  radial  velocity  component  v  is  assumed 
to  be  zero,  so  that  the  radial  momentum  equation  reduces  to 


Figure  2  il'tstraies  the  thrust  producing  components  of  a 
generic  annular  nozzle  The  nozzle  thrust  is  composed  of 
the  thrust  developed  across  the  initial-value  line  FIVl,  the 
thrust  developed  by  the  pressure  acting  on  the  cowl  Fc,  the 
thrust  developed  by  the  pressure  acting  on  the  plug  FP,  and 
the  thrust  acting  on  the  plug  base  Fb  Thus,  the  pressure 
distributions  along  the  cowl  and  the  plug  and  the  base 
pressure  must  be  calculated  to  enable  the  calculation  of  the 
nozzle  performance. 

Swirler  Model 

The  swirler  is  a  set  of  axial  flow  guide  vanes  that  imparts 
the  desired  tangential  velocity  distribution  to  the  air  as  it 
passes  through  the  swirler,  that  is,  from  Section  A  to 
Section  B  in  Fig  1,  Buckley,  Craig,  Davis,  and 
Schwartzkopf1  experimentally  investigated  four  types  of 
tangential  velocity  distributions  at  the  sutrler  exit  They 
are: 


dP/dy  =  pw2/y  (7; 

(c)  The  pressure-velocity  relationship  specified 
by  Bernoulli's  equation  is  assumed  to  be  valid 

dP  =  -  pV  dV  (8) 

This  is  an  approximation  for  the  constant  angle,  the  forced 
vortex,  and  the  Rankine  vortex  swirlers 

(d)  Equations  (7)  and  (8)  are  combined  with  the 
appropriate  swirler  relationship,  Eq  ()).  (2).  (3)  or  (4;,  to 
relate  the  velocity  components  to  a  constant  The  constant 
is  set  by  matching  the  desired  mass  flow  rate  as  described  in 
step  (e)  below.  For  the  free  vortex  swirler,  the  axial 
velocity  component  u  must  be  constant 


Free  vortex: 

* 

n 

n 

y 

(1) 

Constant  angle: 

w-  =  C2u  =  u  tan  a 

(2) 

Forced  vortex: 

u  =  Cjy 

(3) 

Rankine  vortex: 

w  =  (C/y)[l  -  exp(  -  y  Vy j)] 

(4) 

where  the  constants  C,  to  C4  are  specified  to  achieve  the 
desired  tangential  velocity  distribution.  In  the  present 
study,  all  four  swirler  types  were  considered 

For  an  axially  symmetric  combustor  and  nozzle  without 
friction  or  obstructions,  the  tangential  momentum  equation 
reduces  to 

1 5  ■ 


u  =  C5  (free  vortex)  (9j 

For  the  constant  angle  swirler, 

Vvk  =  C6  (constant  angle)  nOi 

where  k  =  sin2a,  and  a  is  the  swirl  angle  in  Eq  (2i  For 
the  forced  vortex  swirler. 

u:  +  2w:  =  C-r  (forced  vortex)  ;|]j 

where  u  is  related  to  y  by  Eq.  (3)  For  the  Rankine  vortex 
swirler.  the  following  integral  relationship  is  obtained 

JCs  v  dv  -  /y  (wvyldy  (12 1 


Dlyu  i/Dt  =  0 


where  CB  is  (he  velocity  magnitude  at  the  inner  radius  and  w 
is  a  function  of  y  as  specified  by  Eq  <4)  The  value  of  CB 
is  determined  iteratively  to  give  the  specified  mass  flow  rate 
at  Station  B 

(e)  The  mass  fiou  rate  is  calculated  using  the 
discharge  coefficient  versus  swirl  results  for  free  vortex 
swirlers  in  Ref  (4). 

(f)  The  mass. averaged  swirl  yw  is  set  by  the 
values  of  the  constants  in  Eqs  (1)  to  (41  Thus,  the  value 
of  Cj  in  Eq  (I)  sets  the  swirl  for  a  free  vortex  swirler;  the 
value  of  C;  in  Eq  (2)  determines  the  swirl  for  a  constant 
angle  swirler;  the  value  of  C3  in  Eq  (3)  determines  the 
swirl  for  a  forced  vortex  sw  irler,  and  the  value  of  C4  in  Eq 
(4)  determines  the  swirl  for  a  Rankine  vortex  swirler. 

(g)  The  mass  flow  rate  at  Station  B  is 
determined  by  numerically  integrating 


m  =  2n  f  puy  dy  (13) 

The  mass  flow  rate  is  a  function  of  the  constant  Cj  in  Eq 

(9)  for  the  free  vortex  swirler.  of  the  constant  C6  in  Eq. 

(10)  for  the  constant  angle  swirler.  of  the  constant  C7  in  Eq. 

(11)  for  the  forced  vortex  swirler.  and  of  the  integration 
limit  Cg  in  Eq.  (12)  for  the  Rankine  vortex  swirler. 
Iteration  for  the  constants  is  required  to  match  the  mass 
flow  rate  in  Step  <e). 

(h)  Equation  (6)  is  numerically  integrated  to 
determine  the  swirl  intensity  yw  Convenient  values  of  yw 
are  obtained  by  iterating  on  the  constants  Cj  to  C*  in  Eqs. 
(1)  to  (4). 

(i)  A  table  of  yw  versus  fractional  mass  flow 
rate  at  Station  B  is  calculated.  As  required  by  Eq  (5).  that 
relationship  between  yw  and  fractional  mass  flow  rate  also 
applies  at  all  other  sections  in  the  flow,  and  specifically  at 
Station  D.  That  relationship  provides  the  needed  initial 
conditions  for  the  transonic  flow  computations  described 
below. 

(j)  The  transonic  flow  calculations  provide  an 
independent  determination  of  the  choked  mass  flow  rate 
through  the  nozzle.  If  the  mass  flow  rate  determined  from 
the  transonic  flow  calculations  does  not  match  the  value 
assumed  in  Step  (e),  then  Steps  (ci  to  (i)  can  be  iterated  to 
convergence. 

To  provide  a  basis  for  comparison  of  different  swirlers,  the 
swirl  number  S  is  employed,  which  is  defined  as 


V;  J  pu:y  dy 


Equation  (14)  is  the  compressible  flow  version  of  the  swirl 
number,  and  can  be  numerically  integrated  to  compute  the 
swirl  number  at  Station  B.  or  a;  any  station  in  the  nozzle 
The  numerator  m  Eq  f!4i  is  closely  related  to  the  mass 
averaged  swirl  yw  in  Eq  (hi.  and  is  constant  at  every 
station  in  the  flowfield  The  disadvantage  of  using  the  swirl 
number  in  this  analysis  is  that  the  denominator  (normalizing 


integral!  in  Eq  (14)  vanes  due  to  changes  in  both  Yj  and 
the  axial  velocity  component  u,  which  changes  due  to  both 
area  change  and  heat  addition  Therefce.  the  swirl  number 
vanes  rather  widely  between  stations  in  the  flowfield  and 
care  must  be  taken  when  using  the  swirl  number  to  compare 
flows  ot  to  correlate  data  to  insure  that  the  values  refer  to 
the  same  station  in  the  flowfield 

Flow-field  Model 

The  accuracy  and  efficiency  of  nozzle  performance 
prediction  techniques  are  highly  dependent  on  tnc  flowfield 
model  chosen  to  represent  the  actual  flowfield  In  the 
present  study,  the  flowfield  model  is  based  on  the  following 
assumptions  O)  steady  axisymmetric  flow,  (2>  tnvivcid 
nonconducting  fluid,  (3)  no  body  forces,  (4)  thermally  and 
caloncally  perfect  gas,  and  (5)  separation  and  base  pressures 
predicted  from  empirical  correlaiions 

The  governing  equations  consist  of  the  continuity  equation, 
the  component  momentum  equations,  the  energy  equation, 
and  the  thermal  and  caloric  equations  of  state  A  detailed 
discussion  of  these  equations  is  presented  in  Ref  (6) 

The  choice  of  a  numerical  method,  or  methods,  to  solve 
these  governing  equations  is  also  crucial  to  ihe  success  of 
the  performance  prediction  procedure  A  subsonic'transomc 
sohRjrm  is  required  to  define  the  mass  flow  rate  and  throat 
thrust.  while  a  supersonic  solution  is  required  to  define  the 
cowl,  plug,  and  base  thrusts  Experience3  has  shown  that 
snbsomc/transonic  flow-fields  can  be  accurately  determined 
fnocn  the  steady  state  solution  to  finite  difference  time 
marching  solutions,  whereas  the  method  of  character istics  is 
the  most  accurate  method  for  determining  supersonic 
nowftelds  in  nozzles.  Therefore,  the  numerical  methods 
.xrJrnnd  for  the  performance  prediction  procedure  in  the 
preweat  study  are  a  finite  difference  time  marching  approach 
totdeftne  the  subsonic/transonic  flow-field,  and  a  method  of 
characteristics  approach  to  define  the  supersonic  flowfield 

SabsoniciTransonic  Flowfield 

lo  the  performance  prediction  methodology  developed  by 
Conley.  Hoffman,  and  Thompson3  for  annular  nozzles 
without  swirl,  the  VNAP  code  (Viscous  Nozzle  Analysis 
Program,)  developed  by  Cline7  was  used  to  calculate  the 
subsonic-transonic  flow-field  In  the  performance  prediction 
methodology  developed  by  Kornblun.  Thompson,  and 
Hoffman4,  the  VNAP  code  was  modified  to  include  the 
effects. of  swirl.  The  resulting  program  wa<  called  SNAP 
(Swirling  Nozzle  Analysis  Program)  Both  the  VNAP  code 
aad  tiie  SNAP  code  are  restricted  to  an  equally  spaced  grid 
«  the  axial  direction  This  is  a  disadvantage  for  nozzles 
having  small  throat  radii  of  curvature,  which  is  generally  the 
e»f*  in  annular  propulsive  nozzles 

In  the  present  investigation,  the  finite  difference  code 
developed  by  Marcum8  was  used  to  calculate  the 
mbtomedran  sonic  flow-field  That  program  solves 
unsteady,  three-dimensional,  mviscid  flowfields  in  super- 
elhptical  nozzles  with  centerbodies  An  option  exists  for 
calculating  unsteady  axisymmetric  mviscid  flow  fields  m 
nozzles  having  circular  cross-sections 

The  gas  dynamic  model  in  Marcum's  analysis  is  based  on 
the  following  assumptions  (1)  unsteady  three-dimensional 
fiou.  (2)  mviscid  nonconducting  fluid  with  no  body  forces, 
and  (3i  a  simple  system  in  thermodynamic  equilibrium  lr 
the  present  analysts,  the  fluid  was  assumed  io  be  g  thermally 
and  caloncally  perfect  gas 
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The  analysis  is  based  on  ihe  unsteady  Euler  equations  The 
finite  difference  algorithm  uses  MacCormack's  explicit  finite 
difference  method'at  interior  points  and  Rentier  s  method10 
at  all  boundary  points  Consequently .  the  method  is  capable 
of  predicting  both  unsteady  and  steady  floss  fields  Steady 
state  floss-fields  are  obtained  by  marching  in  time  until  the 
solution  converges  to  a  steady  state 

The  x>8  physical  space  is  transformed  to  the  iniii 
orthogonal  equally  spaced  computational  space  by  the 
following  algebraic  coordinate  transformation,  which  works 
quite  well  for  a  wide  variety  of  internal  flows 

J  =  arctan  x  (15) 


y  -  Yjt  x  ,0) 

Y 2(x. Ot  -  Yj(x.f)) 


(16) 


a-  =  o 


(17) 


where  Y2(x,8)  is  the  cou-l  radius  ?t  a  given  axial  location  for 
a  given  0  ray,  and  Yj(x,0)  is  the  corresponding  plug  radius. 
The  axial  transformation  permits  the  clustering  of  grid  lines 
in  the  throat  region  of  the  nozzle  For  an  axisymmetric 
nozzle,  Yjlx.B)  =  Yifx)  and  Y|(x,0)  =  Y2(x).  The 
physical  grid  used  in  the  present  study  is  the  81x21  grid 
illustrated  in  Fig  3. 

The  initial  data,  at  time  t  -  0.  for  the  velocity  magnitude, 
the  pressure,  and  the  density  are  obtained  from  a  steady 
quasi-one-dimensional  isentropic  flow  analysis.  The  initial 
values  of  the  flow  angles  are  obtained  by  linear 
interpolation  between  the  centerbody  and  the  wall.  The 
tangential  swirl  distribution  at  the  nozzle  inlet  is  specified  as 
a  function  of  the  fractional  mass  flow  distribution,  yw  = 
yw(rh).  This  distribution  is  determined  from  an  analysis  of 
the  swirler  at  Station  B.  This  swirl  distribution  is  used  to 
compute  the  swirl  yw  at  each  point  across  the  nozzle  inlet  in 
the  following  manner.  After  each  time  step,  the  mass  flow 
rate  is  integrated  across  the  nozzle  inlet  to  determine  the 
fractional  mass  flow  rate  at  each  radial  station.  The 
resulting  fractional  mass  flow  rate  is  used  to  interpolate  for 
the  swirl  yw  at  the  corresponding  computational  grid  pnint 
by  quadratic  interpolation  of  the  swirl  function  yu(rh) 
computed  at  Station  U  The  swirl  velocity  w  is  computed  by 
dividing  the  interpolated  value  of  the  swirl  yw  by  the  value 
of  y  at  the  computational  grid  point 

Steady  state  convergence  is  assumed  when  the  maximum 
percent  change  tn  the  flow  properties,  between  successive 
solution  surfaces,  is  less  than  0.001. 

Points  on  an  initial-value  line  along  which  the  projections  of 
the  Mach  number  in  the  xy  plane  are  supersonic  are  then 
specified  across  the  throat  region.  The  initial-value  points 
are  carefully  chosen  so  that  the  initial-value  line  is  spacelike 
at  all  points  Iteration  on  the  location  of  initial-value  points 
is  sometimes  required  Values  of  the  projected  Mach 
number,  flow  angle,  and  tangential  velocity  along  that  line 
are  determined  by  interpolation  Quadratic  bivariate 
interpolating  polynomials  arc  fit  to  the  nine  closest  grid 
points  to  each  initial-value  line  pom!  In  this  manner,  a 
supersonic  initial-value  line  is  determined,  which  can  be 
used  to  inmate  the  solution  of  the  supersonic  flowficld  by 
the  method  of  characteristics 


Supersonic  Flow  field 

The  supersonic  flowf.eld  ,s  calculated  by  the  method  o| 
characteristics  lor  steady  tw  i-dimensiona!  (low  t nc 
derivation  of  ihe  charac -ensue  and  compat.bilm  equations 
■'  a  «raighttorward.  although  somewhat  lengths,  p?I,Ccduu. 
The  derivation  is  presented  by  Kornblum  and  Thompson ;: 
The  results  are  not  presented  here  In  brief,  three 
characteristics  arc  determined  in  the  xv  plane  the 
streamline  and  the  right-  and  left-running  Mach  bnes 
Three  ctmip.,  .(.nitty  equations  are  determined  along  the 
streamline  an, I  one  compatibility  equation  ,s  determined 
along  each  Mach  line  These  five  conp3tibili(\  equations 
arc  used  to  solve  for  u,  s,  w.  p.  and  p  The  method  is 
general  so  that  cny  consistent  initial  sun!  distribution  can  be 
treated. 


A  direct  marching  method  of  characteristics  procedure  is 
employed  to  solve  the  characteristic  and  conpatibilty 
equations  The  solution  procedure  constructs  the  Mach  line 
network  illustrated  in  Fig  4  The  characteristic  and 
compatibility  equations  are  integrated  numerically  by  the 
second-order  accurate  modified-Euler  predictor-corrector 
method.  The  solution  is  initiated  from  the  supersonic 
initial-value  line,  line  TT‘  determined  from  ihe  transonic 
flow  analysis  The  projected  Mach  number  i  in  the  xy- 
pUne)  along  the  space-like  initial-value  line  must  be 
supersonic  Left-running  Mach  lines  nre  initiated  from  iine 
TT’.  starting  at  point  T,  and  continued  across  ihe  fiowfield 
until  they  intersect  the  cowl,  contour  TE  The  last  such 
Mach  line  from  point  T’  mtersc  :ts  the  cow  l  at  point  T" 
Left-running  Mach  lines  are  then  initiated  from  points  on 
the  centerbody  and  continued  across  the  flow  field  to 
intersect  the  nozzle  wall  This  procedure  is  continued  until 
the  nozzle  exit  l-.p  point,  point  E,  is  determined  Left¬ 
running  Mach  lines  are  then  terminated  at  the  last  right- 
rnnning  Mach  line  emanating  from  Point  E.  This  procedure 
is  continued  until  point  F.  at  the  end  of  the  plug,  is 
determined.  When  Mach  lines  of  the  same  family  coalesce, 
embedded  shock  waves  are  formed  In  the  computer 
program,  when  Mach  lines  coalesce,  one  of  the  Mach  lines 
is  discontinued  This  procedure  works  well  for  weak  shock 
waves. 

The  base  pressure  on  the  face  of  the  truncated  plug  is 
computed  from  the  empirical  correlation  proposed  by 

Johnson1* 


Pt»«=  0  846  Pf/Mf  ■'  (18) 

It  should  be  noted  that  Eq  (I8i  only  fits  measured  data 
when  the  base  area  is  small  It  gives  unrealistically  high 
values  for  highly  truncated  plugs  ansi  for  expansion- 
deflection  nozzles 

The  above  procedure  determines  all  of  the  data  required  to 
evaluate  the  thrust  components  illustrated  in  Fig  2 

R  esuits 

Nozzle  flow  fields  were  computed  for  the  13  cases  defined  in 
Table  1  Figure  5  detines  the  flowtieid  geonetrv  tor  all  13 
cases  The  swirl  strength  was  charas'teri/cd  hs  selected 
values  of  the  mass-averaged  sw  irl  y"w  \  u I u e s  ot  y“u  of  0 
(the  no  swirl  case).  50.  100  and  150  tr.scv  tot  the  Irce 
vortex  swirler,  of  50.  100.  ISO  and  I "5  2r  vc,  (or 
constant  angle  swirler,  of  50.  100,  amt  125  fr  sec  for  the 
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forced  vortex  swirler.  and  of  50  and  100  ftvsec  for  the 
Rankme  vortex  swirler  were  analyzed  Stagnation 
temperature  and  pressure  at  the  nozzle  inlet  were  2500  R 
and  35  psia,  respectively  In  all  cases,  the  specific  heat  ratio 
y  =  1 .40  and  the  gas  constant  R  =  53  35  ( ft-lhl  i  ( Ihn-R  i 
The  values  of  swirl  chosen  represent  the  practical  range  lor 
the  geometry  and  temperature  chosen  here  Higher  values 
of  swirl  result  in  supersonic  Mach  numbers  at  the  suirler 
exit  and  were  therefore  not  considered 

The  compressible  swirl  number  S.  defined  h\  Eq  (14j, 
varies  with  axial  location  in  the  flow  field  Table  I  gives 
approximate  values  of  S  at  Stations  B,  D  and  E  as 
determined  by  numerical  integration  of  Eq  (14)  for  each 
case. 

Figure  6  presents  constant  Mach  number  contours  in  the 
transonic  region  for  Case  I  (no  swirl)  as  computed  by  the 
transonic  flow  analysis  Comparison  with  the  constant 
Mach  number  contours  for  the  other  12  cases  show 
insignificant  differences  Thus.  Fig  6  is  representative  of 
the  plots  for  all  13  cases  There  are  some  minor 
differences,  however,  which  are  discernible  in  the  integrated 
mass  flow  rate  values  and  thrust  values 

Figure  7  compares  the  static  pressure  distributions  on  the 
plug  for  Cases  1,  4,  8,  11,  and  13  There  is  a  pressure 
decrease  on  the  plug  due  to  the  effects  of  swirl  The 
decrease  is  greatest  for  the  free  vortex  swirler  and  least  for 
the  forced  vortex  swirler  The  pressure  distributions  for  the 
intermediate  values  of  suirl  for  each  swirler  design  lie 
between  the  results  presented  in  Fig  7  in  a  very  orderly 
fashion 

Figure  8  presents  the  pressure  distributions  on  the  cow-! 
The  pressure  distributions  on  the  cowl  are  almost 
indistinguishable  for  all  13  swirl  values.  Houever.  there  are 
some  minor  differences  in  the  values  of  cowl.  plug,  and 
base  thrusts. 

Table  2  summarizes  the  mass  flow  rates,  discharge 
coefficients,  and  specific  impulses  for  the  13  cases  studied. 
The  reference  (one-dimensional)  mass  flow  rate  is  21  984 
lbtn/sec.  Mass  flow  rate  values  were  determined  by 
integrating  across  the  initial-value  line  that  is  generated 
from  the  results  of  the  transonic  analysis.  The  zero  swirl 
discharge  coefficient  (0.9853)  reflects  the  two-dimensional 
nature  of  the  flowfield  and  is  a  function  of  the  transonic 
geometry.  The  discharge  coefficients  are  presented  in  Fig 
9  as  a  function  of  mass-averaged  swirl  yw.  The  values 
decrease  with  increasing  swirl  in  a  smooth,  consistent 
manner  and  correlate  very  well  for  all  four  swirler  types 
The  results  indicate  a  high  degree  of  consistency  in  the 
initial-value  data 

Figure  10  presents  the  vacuum  specific  impulse  as  a  function 
of  mass-averaged  swirl  The  results  indicate  that  the 
decrease  in  vacuum  specific  impulse  is  a  function  of  swirler 
design,  and  is  least  for  forced  vortex  swirlers  and  greatest 
for  free  vortex  swirlers  In  fact,  the  decrease  for  the  forced 
vortex  swirlers  is  very  small  This  result  suggests  a 
preference  for  forced  vortex  swirlers  as  far  as  nozzle 
performance  is  concerned,  all  other  factors  being  equal 
The  specific  impulses  of  the  constant  angle  and  Rankme 
vortex  swirlers  he  between  the  values  for  the  forced  and 
free  vortex  swirlers  Notice  that  the  vacuum  specific 
impulse  scale  is  large,  so  that  the  decrease  in  the  worst  case 
is  only  0.35  percent 


In  mosi  cases  of  interest  for  ramjet  or  turbojet  applications, 
the  values  of  the  mass-averaged  swirl  yw  are  no  more  than 
about  100  f!*'scc,  and  the  effect  of  swirl  on  the  nozzle 
periornancc  is  small  and  can  probably  be  neglected  m 
preliminary  design  and  overall  performance  calculations 

The  individual  thrust  contributions  from  the  initial-value 
line,  the  cow],  the  plug,  and  the  base  (illustrated  in  Fig  2> 
are  tabulated  separately  in  Table  3  The  thrusi  values  are 
computed  from  the  method  of  characteristics  (MOC) 
analysis  and  all  values  are  referenced  to  zero  ambient 
pressure 

The  following  observations  regarding  the  effect  of  swirl  in 
annular  nozzles  can  be  made  from  the  resulis  obtained  in 
this  study 

1  The  data  are  extremely  consistent,  giving  a  high  degree 
of  confidence  in  the  trends. 

2  The  effects  of  swirl  on  discharge  coefficient  correlate 
well  with  the  mass-averaged  swirl  yw  for  the  different 
swirlers 

3  The  effects  of  swirl  on  the  vacuum  specific  impulse  are  a 
function  of  swirler  design  The  loss  m  vacuum  specific 
impulse  is  least  for  forced  vortex  swirlers  and  greatest  for 
free  vortex  swirlers. 

4  For  those  applications  for  which  yvT  <  !00  ft*/sec,  the 
effect  of  swirl  on  the  mass  flow  rate,  the  thrust,  and  the 
vacuum  specific  impulse  is  small  and  can  probably  be 
ignored  in  preliminary  design  and  overall  performance 
considerations. 

5  At  higher  values  of  swirl,  the  overall  thrust  is 
significantly  decreased  by  the  effects  of  swirl  In  some 
cases,  part  of  the  thrust  decrease  can  be  recovered  by 
increasing  the  throat  area  to  compensate  for  the  decreased 
mass  flow  rate.  When  geometric  constraints  restrict  the 
nozzle  exit  area,  such  a  design  change  may  not  produce  the 
overall  desired  result,  fn  any  case,  there  is  a  real  loss  in 
specific  impulse  that  cannot  be  recovered. 

6.  W'hen  the  cowl  and  plug  surfaces  are  contoured  to 
optimize  the  nozzle  performance  as,  for  example,  in  Ref. 
(12),  the  optimum  contour  will  be  affected  by  the  presence 
of  swirl  in  the  flowfield. 

Conclusions 

An  analytical  performance  prediction  methodology  for 
predicting  the  effects  of  swirl  in  annular  propulsive  nozzles 
has  been  developed.  Computed  results  for  a  series  of  swirl 
values  for  free  vortex,  constant  angle,  forced  vortex,  and 
Rankme  vortex  swirlers  are  extremely  consistent  and 
indicate  a  high  reliability  in  the  results.  Swirl  decreases  the 
discharge  coefficient,  the  thrust,  and  the  vacuum  specific 
impulse  The  decrease  in  discharge  coefficient  correlates 
with  the  mass  averaged  swirl  yw  for  all  four  types  of 
swirlers  The  decrease  in  vacuum  specific  impulse  is  a 
function  of  the  swirler  design  and  is  least  for  forced  vortex 
swirlers  and  greatest  for  free  vortex  swirlers  For  values  of 
swirl  often  encountered  in  ramjet  and  turbojet  applications, 
the  effect  of  swirl  on  the  nozzle  performance  is  small  and 
can  probably  be  neglected  in  preliminary  design  and  overall 
performance  calculations  Experimental  verification  of 
these  conclusions  is  needed 
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Table  2  Comparison  of  mass  flow  and  specific  impulse. 
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Figure  3.  Typical  subsomc'lranxonic  flowfield  grid 


Figure  4.  Typical  left-running  Mach  line  pattern 


Figure  1.  Geometric  model 


Figure  5  Nozzle  geometric  specifications 


Figure  2  Thru*-:  producing  components  of  an  annular 
nozzle 
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Plug  static  pressure,  psla 


Figure  6.  Transonic  Mach  number  distribution  for  Case  I 


Figure  7,  Pressure  distributions  on  the  plug 


Figure  8  Pressure  distributions  on  the  cow) 


Figure  9  Discharge  coefficients 


Specific  impulse  I$p,  Ibf-sec/lbm 


•  No  swirl 
x  Free  vortex 
o  Constant  angle 
o  Forced  vortex 
o  Rankine  vortex 


Figure  10.  Specific  impulse  comparisons. 
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